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GEJS^VG  is  a  high  energy  laser  propagation  computer 
program  icx  ground- tc-space  applications.  Written  by  Er-  C. 
E.  Hccge  from  the  Air  force  Weapons  Laboratory,  Kirtlanc 
£PB,  it  is  cne  in  a  ianiiy  cf  propagation  cedes  addressing 
this  application-  Specifically ,    GU1SAVG  was  designed  tc 
compute  iriadiance  at  the  target  given  a  model  atmosphere/ 
laser  device  parameters,  and  simple  target  engagement  geon- 
etry.  The  transmitter  induced  effects  or  ream  quality  and 
jitter  are  considered  as  are  the  linear  atmospheric  effects 
cf  scattering/  alscrpticn,  and  turbulence.   A  thermal 
tlocning  model  is  alsc  included.  Adaptive  optics  compensa- 
tion can  he  applied  \ith   consideration  given  to  iscplanatic 
effects. 
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I-  IN5  5CLDC1ION 

GC1SAVG  is  a  simplified  laser  propagation  program.  It  is 
intended  specifically  for  near  vertical  ground-to-space 
applications  utilizing  a  fixed  earth-based  transmitter 
directed  at  a  siigle  target  satellite.   The  primary  purpose 
cf  the  program  is  to  provide  irradiance  and  fluence  en 
target  alcr.g  with  related  propagation  data. 

lie  program  was  written  tv  Charles  E.  Hogge  from  the  Air 
Force  Reapers  Laboratory,  Kirtland  Air  Eorce  Sase, 
Alr>uguer  gue ,  New  Mexico.  GU1SAVG  is  one  of  a  family  cf 
ground- tc-space  propagation  programs.  Other  versiors  irclude 
GOTSIJE  (footprint)  ard  GUTSM1F.  GUTSFP  computes  an  engage- 
lent  envelope  based  or  user  supplied  irradiance  threshold 
levels.  She  propagation  calculation  methods  are  identical  to 
GUTSA^G.  GDXSMIF  is  a  full  wave-optics  program  using  fast 
Fourier  trarsforns  ir  the  beam  propagation  computations. 

lie  lasic  approach  used  in  GOTSAVG  is  to  utilize  modula- 
tion transfer  functions  to  characterize  effects  such  as  team 
guaiity,  jitter,  turbulence  and  to  apply  these  effects  at 
the  aperture  as  a  single  phase  screen-  The  same  approach  is 
used  in  FSf-IV  £Eef.  1].   The -general  modulation  transfer 
function,  (KTF) ,    usee  in  this  way  is  described  in  following 
sections  along  with  the  development  of  each  beam  degrading 
mechanism.   A  linearized  model  is  used  for  thermal  blccming 
which  is  also  effectively  applied  at  the  aperture  as  a  phase 
screen.  In  the  case  cf  blooming,  nowever,  the  phase  variance 
due  tc  blcoming  is  computed,  and  the  Strehl  relation  is  used 
to  detemine  the  relative  irradiance  reduction.   The  effect 
cf  thernal  tlooming  and  the  other  spreading  effects  are 
combined  using  both  I£S  and  multiplicative  methods.  Finally, 
the  average  of  these  two  methods  is  used  to  determine  the 
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total    =_y£t€n    irradiaice    reducticc    xrom    the    diif racticr. 
limited    case. 


1 1 


II-    JHOPAGAIICN    FEATURES 

A.       111    til    APPECACH 

GE15AVG    uses   HTF s   tc   apply    and    taen    to    consolidate    the 
effects    cf    "itteiy    bean    quality,    and    turbulence. 
Essentially,   these    effects    are   replaced    by    a    phase    screen   at 
the    aperture   which    multiplies    the    initial   complex    aperture 
districu ticn  by    a   raiccm   phase    distortion    factor. 

lutcnirsJd.    and    Ytra   £Bet.    2]    have    developed    an    expres- 
sion   icr    the   average    intensity    at    a   point    P    in   the    far-field 
due    tc    perturbations    at    the    aperture.       Ine    intensity    at    P   is 


I(p) 


2ttz 


/ 


M    (f.  ,  z)  exp 


l% 


1 

dp 


/ 


+   I  U(r+bP) U    (r-^p)exp 


i-kirp 


2 

dr 


(2.1) 


uher  e 


p    =     |r2    -    r i 


(2.2) 


and 


r2  + 
r    =    — - -7=r 


(2.3) 


p~  is    the    vector    from    the    z    axis    ox"    symmetry   in   the    far-field 
to    the    tcint  P.       U{r)    is    the    aperture    distribution    and    KJP 
,z)    is    the    M1F    for    tie    disturbance.       c1   and    r2   are    t*o    arbi- 
trary   ccints  in    the    arerture    *here    the    phase   perturbation   is 
ceasuied-    P*  is    the    distance    between   the    points.       z,    in    the 


case  considered,  here,  is  a  constant.      If  the  phase  distur- 
lance  is  a  Landcz    Gdissian  variable  with  a  kncwn  correlation 
functicn,  then  the  Mlf  can  he  expressed  as 

M  (p)  =<exp[i(rj  -r2)]>  (2.4) 

In  terns  of  the  structure  f uncticn, 


M  (p)  =  exp[(-35)D([>(r1-r2)]  <2*5) 


Ihe   seccr.d    integral    in   equation    2.1    represents   the    uircrical- 
ized    aperture    'All.    N  dualizing    this   term    uitn    the    pc^er   in 
the    aperture,     (P0) #    tc   cause    the    flTF    to    be    unity    at    the 
origin   results    in 

Kp)     =     (2^)po/M(|)("p)Ma(p)e"[ikPrfd|  (2.6) 


Kotinc  the  symmetry  of  the  intensity  for  a  given  p  acd 
expressing  M  (P)  as  the  combined  effects  of  jitter,  bean 
Quality,  and  turbulence,  equation  2.6  can  he  rewritten  as 


Kp)  = 


)  =  (^Po/Mj(p)Mt(p)Mb(p)Ma(p)J0(^H)pdl 


(2.7) 


A   Fcurier-Eessel  transform    has   been    used. Note    that    the    fills 
cf    jitter     {&j)t    beam    guaiit y  (M^)  ,    and    turbulence (Mt)    have 
teen    substituted   for    tL^. 

lie  ESJr-IV  manual  ££ef.  3]  contains  the  development 
above  in  mere  detail-  MIFs  for  the  specific  effects  are 
describee    in  fcllowirc    sections. 
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E.       lEESiAI    BLCOflING 

lie    thermal    rlcoiing   model    in   GUTSAVG    is    hased    or.    tie 
following    Jinearized    density    perturbation   equation. 

x' 


Ap_;     *<*>: 


Po 


V 


Y-l    1       f 

Y       poy  ' 


I (x1 ,y) dx* 


(2.8) 


-co 


X    exp  I-    /[<*(£) +o(£)]  d£ 


Here,  a { £)  is  the  aticspheric  absorption  coefficient  at  a 
distarce  £  along  the  learn,  V 0  is  a  constant  transverse  kind 
velocity,  and  p0  is  the  ambient  pressure.  Ihe  exponential 
term  represents  the  tctal  extinction  due  to  scattering  and 
absorptj.cn.   x'is  a  ccnstant  of  integration.  It  is  assuired 
that  the  beam  is  propagated  in  the  positive  z  direction  a 
distarce  £  and  that  the  wind  vector  is  m  the  positive  x 
direction.  Ihe  irtensity  integral  represents  the  heating  cf 
the  atacsphere  as  it  transits  the  beam  £Ref.  4]«  Seme 
assumptions  embodied  in  the  abcve  eguation  are  that  V0<<  c0 
,  the  lccal  sonic  velocity,  so  that  the  process  represented 
cccurs  at  constant  pressure,  and  that  the  kinetics  ci 
absorpticn  and  conversion  tc  heat  are  extremely  fast 
£fief.  5]- 

Ey  applying  the  Gladstone-Iale  relation, the  density 
relation  car  be  expressed  as  a  change  in  the  refractive 
index . 


/    i  x  Ap 

An  =  (n0-l)-^ 
P  o 


(2.9) 
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Eguaticn    2.£  can   be    rewritten    as 

,x' 


An    =    (n0-l)4-  *~  S 
v  o  Y     .  P 


Pu7 

oo 

X      exp|-    Z[aU)+aU)]  d£ 


(2.  10) 


Ihe  chance  in  the  wavefrcnt  phase  of  a  team  due  to  the 
refractive  index  charge  wner  the  beam  is  propagated  a 
distance  i  ,    is  given  iy 


A4> 


2_TT  /* 


i 


And£ 


(2.  11) 


£ubs titu tin^  equation    2.10    intc    2-11    and    changing    the    limits 
cf    integration    tc   reflect   the    ground-to-space   propa^atrcn 

path    results  in 

2tt 


x 


A<J>(x,y)     = 


X 


n„-l    Y-l  f 
Po  Y    7_, 


X 


I (x' , y) dx* 

h 


atm 
a(h)exp| -sec6   /   [a  (h)  +a  (h)]  dh 

ht 


\2.  1z) 


dh 


V0cos^    +    coh 


h    ,       is    the    extent    cf    tie   atmosphere,    which   is   about    atcut 

atm  K 

30  km,  aid  h,  is  the  height  of  the  laser  transmitter.  Alsc, 
the  wind  term  has  beer  expanded  to  include  tne  relative  wind 
velocity  due  to  slewing.  V  0  is  assumed  to  be  parallel  and 
opposite  in  directicr  of  that  cf  the  target  mction.  E,   is  the 
angle  cf  incidence  oi  the  wind  to  the  beam  so  that  V0ccs{Q 
represents  the  transverse  wind,  coh  is  the  effective  wind 
generated  ly   slewing.   co  is  the  angular  slew  rate. 

Ihe  first  pcrtion  of  eguation  2.12  is  independent  cf 
path  while  the  second  part  is  not,  assuming  I(x,y)  dees  net 
change  alcne  th.e  propagation  path.  This  assumption  .is  valid 
enly  fcr  a  very  small  amount  of  blooming-   The  approach  used 
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ir.    GUTSAVG    is    to    determine    the    phase    distortion    1  ue    to 
thermal   blooming    by    first    evaluating    the    path    invariant    part 
of   equation    2.12    .    This   is    accomplished    by    constructing    a 
phase    screen   at    the    aperture    and    then    removing   the    best    fit 
tilt,    focus   curvature,    ani    mean    phase.    ZerniKe    polynomials 
are    used    to    model    these    aberrations.    The    result    is    the    resi- 
dual   phase    due    to   thermal   blooming    alone.       The    variance    of 
the    phase    is   then   compute:!. 

The    path   dependent    term    is    evaluated    within    the    angle 
interval   loop   of    the    program    and    is    applied    to   the    previ- 
ously   computed    phase    variance    during    each    path    iteration. 
The    path    iteration    process    is    diagrammed    in    the    engagement 
geometry   section.    Also,    see    Figures    2.1    and    2.2    for    a    flow 
diagram    of    the    general    treatment    of    thermal    blooming    ir    the 
program. 

Once   the   total    phase   variance   has    ueen    determined,    the 
Strehl  relation    is    used   to    compute    the    intensity    degradation 
due    to   thermal    blooming. 


i  =    exp(-a2)  (2-  13) 


The    result    of   equation    2.13    is   a    relative    intensity     (Irel) 
ratio.    Io   is  the   ideal   on-axis    irradiance    with   no    phase 
distortion.      The    Strehl   relation    above    is    thought    to    be    too 
severe   a   model    for    Irel   below    0.3    [Ref.    6].       For    that 
reason,    if    a2     is   less   than     1.2    ,    the    Irel    will    oe   computed 
using    polynomial    curve   fits    developed    from    GUTSMTF    results. 
GUTSMTF    is    a  full   wave  optics   code    utilizing    fast    Fourier 
transforms.      For   a    description    of    the    curve    fit    method 
above,    see    the    subroutine   BLOOM    explanation.       [Ref.    7] 

Combining    the   thermal  blooming   effect    with   the    other 
effects,    such   as   turbulence    and    jitter,    is    accomplished   by 
averaging   the   results   of    two   different    approaches.    The    first 
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approach  is  the  RSS  (xoct  sum  squared)  method.  Ihis  method 
ex  ccmriiing  tne  Irel  due  to  the  effects    ox  thermal  tlccming 
with  the  Irel  due  tc  fitter,  team  dua lxt y ,  and  turbulence  xs 
acccHj:  J  isied  as  follcws 


-l 


Irel 


rss 


=   1+ 


Irel 


-1 


tb 


r  i 


Irel 


-  1 


o 


(2. Mi) 


result  ana  Irel   xs  tie 
tb  o 

result  oi    the  other  eifects.    Ihe  second  approach  is  a  multi- 
plicative approach  arc  is  simply 


Irel 


m 


=  Irel   x  ireL, 
o        tb 


(2-15) 


Ihe  tv«c  confined  Irels  chtained  ty  these  methods  axe  then 
averaged  tc  give  the  total  intensity  ratio  due  to  ail  the 
attenuating  or  distorting  propagation  effects. 


-  =  irel.. 
I  o        tot 


Irel   +  Irel 

m        rss 


(2.  16) 


I   is  the  ideal  diffraction  iiiited  on-axis  xEtensity.   Ihe 
tases  fcr  the  above  averaging  process  is  empirxcal  in  nature 
and  is  ar  attempt  to  adjust  the   results  obtained  ty  the  E£S 
method  alcne.   The  results  produced  by  ESS  were  thought  tc 
te  tec  cj:  tinistic.  ihe    multiplicative  method,  a  mere  pessi- 
mistic approach,  was  therefore  included.  Ihe  ultimate  Irels 
cbtaired  are  very  close  to  these  obtained  ty  GUISMTF,  the 
full  wave  cede.   [fief.  8] 

Ihere  are  sone  linitaticns  tc  the  thermal  blooming  model 
used  in  GG1SAVG  in  addition  tc  the  assumptions  already 
mentioned,  first,  the  wind  (V0)  is  applied  as  a  corstant 


M 


everykteie  in  the  atncsphere.  At  nigher  altitudes,  this  is 
ret  a  najcr  consideration  due  to  the    higher  relative 
velocity  cf  the  team.   At  lew  to  medium  altitudes  this  cculd 
affect  tiernal  blooming  to  an  extent  to  warrant  the  addition 
cf  a  wind  profile  as  a  functioE  of  altitude.   This  cculd  te 
done  *itl  little  effort  if  the  data  is  available.   Also,  the 
program  assumes  a  wind  parallel  and  opposite  to  the  direc- 
tion ci    the  target  mction.  This  precludes  the  case  cf 
slewirg  kith  the  wine  and  the  creation  cf  null  spots. 
Transcric  blooming,  which  violates  one  cf  the  origiral 
assumptions  of  the  ecuation  used,  is  not  considered 
£Eef.  9]-   for  very  lew  altitude  satellites,  the  high  sleu 
rates  gererated  fcoulc  result  in  a  supersonic  relative  wire 
across  tie  team.  In  this  case,  the  constant  pressure  assump- 
tion  is  irvalid  [Ref.  10]-   Kinetic  cooling  and  molecular 
breakdown  are  also  net  addressed  in  GUTSAVG- 

C.   SCATTERING  AKD  AESCRPTICN 

GC3SAVG  does  not  contain  any  type  of  atmospheric  rcccel 
kith  respect  to  scattering  and  absorption  data.   Extinction 
coefficients  must  be  entered  by  the  user  in  the  appropriate 
subroutines,  ALPS  and  ALfA,  or  the  program  can  be  modified 
to  accept  a  separate  data  file.   Once  inserted  into  the 
program,  scattering  aid  absorption  coefficients  are  treated 
kithcut  distinction  between  aerosol  and  molecular  mecha- 
nisms. Therefore,  the  coefficients  used  must  represent  the 
total  effect  of  scattering  or  absorption. 

a(h)  =   a(h)mol  ♦  a(h)aer  «*•"> 
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cr(h)    =      a(h)mQl    +    CT(h)__  (2.18) 


aer 


"Transmission   due    tc    scattering    and   absorption    are 
computed   indenticall}'   and   given    ty 

Ta  =  exp  -sece  ja(h)dh)  'U- 1S) 

ht 
where  6  is  the  zenith  angle,  ard  the  integration  limits  are 
the  altitude  of  the  transmitter  and  tne  vertical  extent  ot 
the  a tuos there .   The  extent  oi  the  atmosphere  in  the  program 
is  defined  as  30  km.   Figure  2.3  snows  tne  prograz  applica- 
tion cr"  scattering  aid  absorption. 

GllSflvG  allows  tie  user  tc  describe  the  neam  at  the 
aperture  in  terms  of  the  electromagnetic  field  amplitude 
distribution  and  total  beam  quality.  Ihe  intitai  amplitude 
field  used  ty  the  prcgram  is  Gaussian  in  sha^e  with  the  user 
specifying  c  waist  diameter.  The  waist  diameter  is  defined 
ty  the  1/e2  point  on  the  distribution.   Truncation  cf  the 
this  Gaussian  field  will  of  course  depend  en  the  aperture 
diametei  and  the  size  of  the  central  obscuration.  If  the 
waist  diameter  is  made  large  cempared  to  the  aperture,  a 
jiore  uniforn  distribution  results. 

Beam  guality  at  the  aperture  exit  may  be  specified  by 
two  different  parameters.  One  cf  tnese  is  the  'times 
diffraction  limited  number' ,    N.   N  has  been  used  in  a 
general  way  to  mean  an  increase  in  far-field  spot  size  cr  as 
a  ' pewer-in-the-bucket '  ratio.  In  GUTAVG,  N  is  a  total  team 
guality  ten.   The  second  beam  guality  parameter  is  a  nendi- 
nensicnal  term  representing  the  EMS  phase  distortion  at  tie 
laser  wavelength  at  tie  exit  aperture,  ^rms  - 


1S 


3he  paraaeters  axe  related  to  each  ctaex  and  to  the  inten- 
sity degradation  by 


I 
Io 


N 


,'2tt5 
exp  - j    rms 

A 


(2.20) 


1c    apply  beam  guality  tc  the  propagation  problem,  ai  E1F 
array  is  constructed  representing  a  phase  screen  at  the 
aperture.  Ihe  phase  is  assumed  tc  be   a  Gaussian  random  vari- 
able with  a  zero  mean  value.-  Per  the  axi-symmeti ic  lean, 
considered  in  GU1AVG,  the  Ali    is  [fief-  11] 


Mfa(p)  =  exp(-k 


°  -cA(P: 


(2.21) 


C      is    the    autocorrelation   function    of    the    phase    and    is 
defined    as 


v^  = 


o    exp 


-(if 


(2.22) 


2^5 


where 


is  the   phase   variance,     ( 


rms      2  , 

-r )  ,    and    L    is    the 


phase    correlation   length   [fief,    12J.       The    neam   quality    HIP 
array    is   combined    with    the    MTF    arrays    due    to    ether    propaga- 
tion   effects  to    detemine   the    conplete    system    MTF    and, 
hence,    tie    total   irradiance    degradation.       Based    on    the 
provided    input    parameter,    Figure    2.4    shows    the   general 
treatnent    of  beam    guality   within    the   program. 


E-       ICfiEOIFNCE 

GtSS^VG    uses    the    Eufnagel    irodel  ffief.     13]   for    C2    as    an 
indicator    of   the   optical   effects    of    turbulence   along    the 
propagation   path.    C2    is   the    refractive   index   structure 
constant   and  represents    the    refractive   index    in    the 


atmosphere  as   a  furcticn  of  turbulence  induced  density 
fluctuations.  The  model  is  an  empirically  derived  vertical 

profile  ci  C2  . 

n 

Iriec  [Bef-  14]  las    developed  a  parameter  which  is 
directly  related  to  tie  behavicr  cf  a  coherent  beaa  ir  a 
turbulent  medium.  This  term  is  called  the  effective  coher- 
ence ciaaeter,  r  .  Ir  the  case  of  a  laser  trarsmiter,  c0 

o 

represents  a  physical  limit  to  the  transmitter  diameter  oi  a 
rear  diffraction-limited  beam,  ror  a  transmitter  diaaeter,  D 
,    larger    than   r      ,    degradation    cf    the    ceam    b_y    turbulence 


viill    cccir.    If    D    is    saaller    than    r 


tnen    near    dirfracticc- 


limited    propagation    >»ill   be    achieved.    Beam    wander    oi    puit 
tilt    cccurs   for    transnitter    diameters    approximat ely    egual    to 
r    .       Jura    [Sefo     t5  ]    has    defined    a    somewnat    different    but 
related    tern  that   car  be   thought    of   as   a    lateral    cohererce 
lengtl,    p0.      These    t*c    guan titles   are    given    by 


r  (.    = 


l^l  k2sece   /  C^ 


(2.23) 


ana 


Po    = 


1.45    k'sec 


'U 


(h)    dh 


-3A 


(2-24) 


so    that 


Po    = 


-       ro 


2.1 


(2-25) 


inhere    x    = 


2rr 

A 


is    the   zenith    angle,    and    tie    limits    cf 
,    the   vertical   ext 
and    It    #tte    height    of    the    transmitter. 


integration   are    b         ,    the   vertical   extent   cf    the    atmosphere, 

atm 


Po       may    be   specified    by    the    user    as    a    program    input    or 
the    program   may    be    allowed    to    compute    it    based   oi     equation 
2. 24    .    The    user    may    also   input   Po     indirectly    by    specifying 
the    'seeing   conditions'    ,    a    quantity    ased    by    astronomers    to 
describe   angular   spread    of    a   stellar    point    source     (see    input 
section   and   Figure    2.5).       I  he    progran    uses    P0    to    compute    the 
atmospheric   MIF. 

The    MTF   of    the    turbulent    atmosphere   is    determined    by 
developing    the    structure   function      of    the    turbulence.       This 
development   is    demonstrated    Dy    Yura    * Ref .     16].      The 
resulting    MTF    is   given   by 


M    (p)    =    exp 


P_ 

Po 


(2.  26) 


This    MTF   effectively    applies   the    the    effect    of   tirbulence 
along    the    propagation   path    as    a    phase    screen   at    the 
aperture. 

F.       JITTER 

Beam    jitter    is    a    user    input    and    is    specified    by    the 
variance    of   the    angular    excursions    of    the    beam.    J  sing    the 
2-sigraa-p   definition, 


2  a        =  V2(az    +a2) 
P  x        y 


(2.  27) 


where   a      and   av    are    the   axial    variances    of    the   jittered    beam 


x    — -      y 
center   in   the   far- field 


tf      and   aw    are   random    variables   with 

a  y 


Gaussian    distribution  and    in   the    symmetric    case    ,    as   consid- 
ered   by    GUTSAVG,    0^=    a      .    The    resultant    intensity    distribu- 
tion   due   to    jitter    will   also   be   a   Gaussian    distribution   with 
{2a   )     representing    the  spot    radius    defined    at    the    1/e2 
point.    In    other    words,    86.5%    of    the    beam    energy    will   reside 


within    the    radius    2a. 


[Ref.    17] 
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Jitter  can  be  shewn  to  re  a  wavefront  tilt  at  the  aper- 
ture. Using  this  approach,  a  phase  screen  at  the  aperture 
can  te  used  to  characterize  the  effects  of  jitter  and  a 
jitter  MTJ  developed.  That  MTF  is  given  by 


M  .  (  p)  =  exp 


-k2p2(2o  )  2 


(2.28) 


figure  2.6  is  a  general  flow  diagram  for  the  treatment  cf 
jitter  within  the  program. 

G-   AIAE1IVE  OPTICS 

lie  ability  to  apply  adaptive  optics  corrections  to  tie 
propagation  problem  has  bee n 'included  in  GUTSAVG.  The  user 
has  several  options  kith  respect  to  the  type  and  degree  or 
compensation  desired.   The  following  general  discussion  and 
figure  2.8  provides  the  needed  insight  to  the  effects  of 
selecting  the  adaptive  optics  options. 

The  najcr  compensation  mode  provided  by  the  prcgran  is 
invoked  ty  selecting  full  zcnal  adaptive  optics  with  consid- 
eraticE  given  to  isoplanatic  effects.   When  selected,  this 
model  results  in  the  correction  of  beam  degradation  due  to 
turbulence.  This  is  accomplished  by  correcting  turbulence 
induced  tilt  and  then  adjusting  Po  so  as  to  produce  a 
predeternined  level  cf  adaptive  optics  performance.  This 
predeteriined  performance  is  as  measured  by  the  Strehl  ratio 
given  a  residual  phase  variance  determined  by  the  adaptive 
optics  senscr  phase.   Parameters  determining  the  phase  error 
are  the  response  bandwidth  of  the  adaptive  optics  system, 
the  number  cf  system  actuaters,  the  reflected  radiant  inten- 
sity of  the  target,  and  the  target-to-sensor  transmission. 
Ihe  resultant  .Po  found  in  this  manner  is  then  used  to 
compute  the  atmospheric  MTF .  When  this  'corrected'  Mlf  is 
used  to  compute  the  far-field  intensity,  the  result  will 
represent  an  adaptive  cptics  corrected  value. 
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fcithcut  invoking  full  zcnal  adaptive  optics,  the  user 
flay  apply  a  tilt-only  correction  for  turbulence.  In  this 
case  seme  cz   all  of  tie  tilt  due  to  turnulence  may  fc€ 
removed  hefcre  computing  the  atmospneric  i-ilr.  The  degree  cf 
tilt  compensation  is  specified  by  the  user. 

As  mentioned  above,  isoplanatic  effects  are  included  in 
the  adaptive  optics  calculations.   This  is  also  an  cpticn. 
however,  and  isoplanatic  calculations  may  he  inhibited  ry 
the  user.  The  effect  considered  is  the  limitation  cf  the 
adaptive  optics  systen  giver  an  isoplanatic  angle  smaller 
than  tie  target  lead  angle.   Fried  (Bef.  18]  provides  a 
discussion  cf  iscpiaratism  and  development  cf  the  iscpla- 
ratic  arcle. 

Although  adaptive  optics  compensation  ror  thermal 
hiocnirg  is  not  modeled  in  a  strict  sense,  a  thermal 
blooming  correction  factor  can  he  applied.  This  factor  is 
simplv  a  fractional  constant  that   multiplies  the  thermal 
blooming  phase  variance  before  the  Strehl  relation  is  used 
to  cenpute  the  intensity  degradation. 

Eefer  to  the  input  definitions  and  the  subroutines 
involved  with  adaptive  optics.   In  addition,  refer  to  figure 
2-7  fci  ncxe  explanation. 
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Figure  2. 1 


GC1SAVG    1-hexmal    Blooming    Algoritha 
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Figure  2-2   Tlermal  Elccming  Algorithm  (coct) 
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Figure    2.3        GU1SAVG   Scattering    and   Absorption   Algorithm 
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Figure  2-4    GD1SAVG  Eeam  Quality  Algorithm. 
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figure    2-5        G01SA¥G    lurfculence   Algorithm. 
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figure  2.€   GDTSAVG  Jitter  Algorithm, 
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figure  2-7    GL1SA1G    Adajtive  optics  Algorithm 
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Figure  2.8   GOTSA\G  Adaptive  Optics  Algoritha  (cent) 
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Figure  2.9    GOTSflDG  Adaptive  Optics  Algorithm  (cent) 
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Figure  2.10    GOTSATiG  Adaptive  Optics  Algorithm  (cent) 
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Figure  2-11    GUTSAiG  Adaptive  Optics  Algorithm  (cent) 
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Figure  2-12    GUTSiVG  Adaptive  Optics  Algorithm  (cent) 

36 


III.     JPBOGBJJ    AND    SOEffiOGBAfl    DISCfilPlIOM 

A.       IBCGFAfi    INP01S 

Ihe    following   are   tie   ir.puts   necessary    to    utilize 
GUTSAvG.    Default    values    are    indicated    for    parameters    not 
ansclutely    required    fcr    program   operation.       At   NPS    these 
input    tarameters  are    entered   via    an  input   file.       A    copy   of 
this    file    i£   provided    in   Appendix    A.       (*)     indicates    a    ncndi- 
aensicral    parameter. 

•  EIA  (meters) 

Ihe    ciaaeter   of    the    transmitter   aperture. 

•  IIACES    (meters) 

lie    ciaaeter  of    the   central    obscuration   of    the    trans- 
mitter . 

•  EIAMSZ    (meters) 

Ihe  Gaussian  waist  diameter  of  the  amplitude  distriru- 
ticn  at  the  aperture,  (measured  at  the  1/e2  point) 

•  RAvE    (meters) 
Ihe  laser  wavelength. 

•  EICIAL  (Hatts) 

Ihe  total  power  at  the  aperture. 
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Ihis  is  the  cft€i  used  "tines  diffraction  limited 
cuiter"  and  represents  total  team  quality.  As  used  in 
GESS,  it  is  related  to  the  R'AS    phase  distortion  at  the 
ajerture  by  the  £trehl  approximation 


(TDFLMT) 


2=  exp 


2tt5 

rms 


(3.1) 


(ILflMl)2    ,     therefore,    is    equivalent    to    the  c atic    of    the 
cr-aiis    diffraction   limited    intensity    to    the    cn-axis 
intersity   resulting   from    the   near-field    phase    distci- 
ticn,  6 


rms 


(TDFLMT) 


I 
Io 


{3.2) 


Per    a    discussion    cr    the    limitation    cf    3.1,    see 
£Sef.    IS], 

•    fia?EEC    (*) 

Ibis  term  is  the  EMS  phase  distortion  at  the  laser  aper- 
ture, rcndinensicralized  hy  the  wavelength. 


WAVEBQ  = 


rms 


[3.3) 


•  SCAIEQ  (meters)     default  =  DIA/5 

Ihis  term  is  the  transverse  phase  correlation  lerctt  at 
tt«  aperture. 
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•  THSEE    (arcsec)  default    =    f(RHOO) 

This    is    a   qualitative    term    used    by    astronomers    to 
describe    'seeing    conditions'    in    the    visible    range.     If    a 
point   source    is    viewed    from    the    earth,    it    may    not    appear 
as   a    point    source   but    as   a    'smear'    or    spot.    The    angular 
spread    of   this    spot    is    the    parameter   THSEE.    If    RHOO    is 
not    specified   as    an    input   and    THSEE    is,    the    program    will 
use    THSEE    to    compute    EHOO.    Other    than    this   case,    THSEE 
is   not    used    but    will   be    computed    as    a    function    of    RHOO 
and    included    as    an  output. 

•  HGBND    (neters) 

The    height    of   the   ground    at    tne    transmitter    position 
above    MSL. 


•  HTRANS    (meters) 

The    height    of    the   transmitter    above. 

•  HSAT       (meters) 

The  orbital  altitude  of  the  satellite  above  1 SL  (at 
zenith)  . 

•  THETHX  (degrees) 

THETMX  is  the  angle  measured  from  zenith  below  which  the 
laser  will  not  transmit  for  a  zenith  pass.  For  offset 
flight  paths,  it  should  be  noted  that  the  transmitter 
will  point  below  this  value.  For  an  illustration  of  the 
engagement  envelope,  see  figure  3.1 
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•  ICIF   (meters) 

Per  a  target  that  dees  ret  pass  direct  overhead,  this 
input  specifies  tie  amount  that  the  target  ground  ttack 
is  offset  from  the  ground  track  of  the  overhead  case.  It 
is  tie  distance,  as  measured  from  the  transmitter,  to  a 
perpendicular  intersection  en  the  ground  track. 

•  EHCO    (meters)     default  =  f(IriSEZ)  or  SHCI^E 

Ifcis  is  the  turbulence  coherence  length  as  defined  ty  H. 
1 .  i'ura.   £Eef-  2C]-   {see  turnulence  section  fcr  mere 
discussion  aid  reierencts) 

•  VO     (m/sec) 

Ire    atmospheric    fcind.    The    direction    of    the    wl  cd    is 
parallel  and  opposite   to    the    direction    of    target    ncticn, 
i^cte,    in   the   present   program,    VO    is    a    constant    indepen- 
dent  of   altitude. 

•  ACEICfl  (*)     0-0  to  1-0 

Ihis  parameter  allows  the  user  to  correct  for  thermal 
ilccming  as  if  by  adaptive  optics.  I  he  value  entered  nay 
range  from  0.0  tc  1.0.  If  0-0  is  used,complet e  thermal 
rlccning  compensation  will  occur.   Conversely,  if  1.0  is 
entered,  no  compensation  will  he  applied-   The  variance 
ci    the  phase  distortion  is  multiplied  by  this  correction 
factcr  fefore  the  Streh 1  relation  is  used  to  compute  the 
relative  intensity  reduction  due  to  thermal  blooming. 

-   =  exp(-o2  *  AOBLOM)  (3-4) 
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•  AVGSPT  (*)    0..0  to  1  .0 

AVGSIT  allows  the  user  latitude  in  defining  the  far- 
field  spot  size  tc  otner  than  that  indicated  by  GUTSAVG 
analysis-  A  value  of  1.0  would  result  in  the  peak  irra- 
diance  according  to  the  program  analysis.  An  entered 
value  of  .5,  for  example,  would  result  in  a  peak  irradi- 
ance  50?  less  than  the  program  would  otherwise  indicate. 
AVGSIT,  then,  is  an  adjustment  factor  that  allows  the 
user  tc  account  for  effects  not  addressed  in  the  GUISAVG 
propagation  calculations. 

•  SIGJI1  (radians) 

SIGJIT  is  the  2o    variance  for  pointing  and  tracking 
jitter. 

<•  A£AP   {*)    0-  C  to  1.  0 

This  term  is  a  correction  factor  for  tilt  due  to  turbu- 
lence.  It  represents  the  residual  tilt  after  AO  conpen- 
saticn.   ADAP  may  he  varied  from  0.0  to  1.0-  A  value  cf 
1.0  would  result  in  no  tilt  due  to  turbulence  being 
removed  while  a  value  of  0-0  would  result  in  total  tilt 
due  to  turbulence  'compe  nsation.  If  full  AO  is  selected 
tie  program  will  set  ADAP  egual  to  0.0. 

•  NF1AGA  (*)    0  cr  1 

NPIAGA    is  a    selection  indicator   for   full    zonal    adaptive 
cttics.      Enter    a    0   if   AC   is    not    desired,    enter    a    1    if    AO 
is    desired.     If    AC   is   used,    XJT,BWIDTH,    and   NA    must    be 
specified. 
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•  hCISC       [*)         0  cr  1 

isCISC  is  a  selection  indicator  for  isoplanatic  Cclculc- 
ticrs.   Enter  a  1  to  inhibit  isoplanatic  calculatic is ; 
enter  a  0  if  iscplanatic  calculations  are  desired. 

•  XJ3     (Watts/sterad) 

XJ1    is    the    target   radiart    intensity.    It   is    one    of    the 
factors    used   to    determine   how   noise   aifects    the    res^orse 
cf    the    P.O   system. 

«    EEIEIE    (Hz) 

EJiHIH   is    the   barcwidth.    of    the    adaptive    optics    system. 

•  UA      (*) 

H&    j.s  the  nunber  cf  AO  actuaters  used  tc  perform  phase 
ad jtstnent. 

•  AES1C2  (*) 

This  is  tat  percent  transmission  at  zenith  at  the 
sensing  wavelength  of  the  AO  system.  This  term  is  used 
ic  the  deter nination  of  relative  noise  at  the  AO  serscr. 

•  M1      (*) 

Ihe  rumrer  of  iteration  steps  for  the  argle  lccp  cf  tie 
prcgrajj  from  THE1KX  to  0- 
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•  N2      (*) 

lie  numier  of  altitude  intervals  for  absorpti  or.  and 
scattering  deter miration. 

«  M3      (*) 

Ifce    rumrer    of    altitude    intervals   for    the    turt*  ulence 
calculations. 

•  N4      {*) 

lie  numier  of.   iteration  intervals  for  the  MTF  Calcula- 
tions. 

•  M5      (*) 

lhe  rumrer  of  intervals  for  the  slant  path  update  cf 
thernal  blooming- 


E.   ISGAGEHENT  GEOHElfY 

lhe  target  engagement  -geometry  is  that  of  a  earth-rased 
transmitter  and  a  target  satellite  at  a  given  orbit.  Nc 
attempt  is  cade  to  define  an  engagement  envelope  based  en 
threshold  irradiance  cr  fluence1.   The  input  parameters 
defining  the  engagement  window  in  GUISAVG  are  THEI21X,  ICFE, 
and  HSA1.  lhe  general  geometry  of  this  window  is  shewn  in 
Eigure  3-1  .  Only  half  of  the  total  transit  window  is 
addressed  in  the  program  calculations;  the  results  are  the 


1    lhe  GUISFP  (fectprint)  version  of  gtts  was  writter 
to  dc  this.   Except  fcr  this  feature,  the  propagation  calcu- 
lators are  identical  to  GUISAVG. 
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same  fcr  either  half.   The  program  output  reflects  this 
half-window  evaluation  except  for  parameters  such  as  total 
fluence  and  shot  time  which  are  simply  dourle  the  ccsputed 
values- 
Figure  3.2  defines  the  earth  center  angle  (ECANG) .  FCANG 
is  a  function  of  the  user  input  IHEIMX-  Most  of  the  geome- 
tric calculations  are  referenced  to  earth  center.  Therefore, 
this  argle  is  used  fcr  computing  such  positioning  data  as 
the  angle  irterval  at  which  the  irradrance  will  ue  evaluated 
(see  Fi  gure  3.  3)  . 

Cffset  flight  paths  require  a  coordinate  translation  as 
snown  in  Figure  3.4  .  Position  ana  velocity  relative  to  the 
transcitter  are  computed  as  in  Figure  3.5  .  It  should  re 
noted  that  if  the  flight  path  is  offset,  the  zenith  angle 
will  exceed  IHEIMX  fcr  part  or  all  of  the  window.  This  is 
tecause  the  windcw  is  defined  in  the  x-z  plane  only. 
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H   -  Zenith  altitude  of  target 
s  ^ 


m 


-  Maximum  zenith  angle  for  engagement  (zenith  path) 


L   -  Offset  distance 
o 


Zenith 
Flight  Path  i\ 


Offset 

Flight  Path 


Engagement 
Window 


x 


Figure  3.1    General  Engagement  Geometry 
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X 


=  cos    —   sin  6   +0   -  7r 
R       m    m    2 
s 


Target  satellite 

Laser  transmitter 

Earth  center  angle 

Maximum  zenith  angle  for 

engagement 

Earth  center-to-target  distance 

Earth  center-to-transmitter 
distance 


figure  3.2    Earth  Center  Angle  (ECANG) 
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Irradiance  is  computed  at 
interval  midpoints 


The  number  of  intervals  is 
a  user  input  variable 


Figure  3.3    Angle  Intervals  for  Irradiance  Evaluatici, 
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Offset  is  from  T  to  C 


K  =  Rsin(7/2) 

Ly  =  cos ( y/2 ) • 2K 

Lz  =  sin (y/2) • 2K 


figure  3.4   Coordinate  Translation  due  to  Flight  Path  Cttset 
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Velocity  components 


V  =  -V  cos  (9  ) 
x      s      e 

V  =  V  sin (9  ) cos ( y) 
y     s      e      ' 

V  =  V  sin(9  ) sin(y) 

Cj  O  t^ 


Position 

X  =  R  sin( 9  ) 
s      e 

Y  =  Lycos  (y)  +  (Z  o-R-Lz) 

•  sin (y) 

Z  =  -Ly • sin (y) + (Z o-R-Lz 

•  cos ( y) 
where   Z0=  R  cos(9  ) 


Range  from  transmitter (T)  to 
satellite (S) 


R^   =   Z  +  X  +  Y 
ts 


(For  Ly  and  Lz  definition,  see  Figure  3.4) 


figure  3.5    Velocity  and  Position  Relative  to  Transmit ter. 

4S 


C.       HUM    EECGBAfl    F10K    DIAGBAM 

lie    nail  program    iicstly    consists   of    geometry    calcula- 
tions   and    decision    flew    points.    Ihe   decision    foiu ts    allcw 
branciirc    tc  adaptive    optics    and    isopianatic    subroutines. 
Most    cf    the    propagation   calculations    are    done    within    the 
subroutines.      Only    tve   major   iteration    loops   reside    in    the 
aain    jrccran,    the   ancle    interval    loop    and    the   combined    Elf 
integraticr    loop.      Tie    following    is   a    general   flow    dia^rac 
for    tie    nair  program.    Decision    variables    are    user    lifuts    cr 
program    defined    jiaraneters. 
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figure  3-6   flain  Program  Flow  Diagram. 


o 


FAPTLD 

COMPUTE    THE 

FAR-FIELD    IRRAOIANCE 

ANO    THE 

CORRESPONDING 

APERTURE    MTF 


NO 


WAVEBQ-0 


YES 


BOIREL 


COMPUTE    TDFLMT    USING 
THE   SPECIFIED  WAVEBQ 


DETWAV 


DETERMINE  THE  UAVEBQ 
THAT  WILL  RESULT  IN 
THE   SPECIFIED   TDFLMT 


8M 


CONSTRUCT  THE  BEAM 

QUALITY  MTF  ARRAY  AT 

THE  APERTURE 


JIT 


CONSTRUCT  THE  JITTER 

MTF  ARRAY  BASED  ON 
THE  SPECIFIEO  JITTER 


Figure  3.7    Bain  Program  Flow  Diagram  (cout) 
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figure  3-8    Bain  Program  Flow  Diagram  (cont) 


c  z 


SCINT 


COMPUTE  THE  LOG 
AMPIL1TUOE 

SCINTILLATION  FOR 
ZENITH  TRANSMISSION 

THIS  IS  USEO  LATER  FOR 
AO  CALCULATIONS 


PRTPHS 


COMPUTE  THE  THERMAL 

BLOOMING  PHASE  CHANGE 

INDEPENDENT  Of  PATH 


RESPHS 


COMPUTE    RESIDUAl     PHASE 

OUF    TO    THERMAL 

BLOOMING   BY    REMOVING 

MEAN.    TILT,    ANO 

FOCUS   CURVATURES 


COMPUTE    THERMAL 
BLOOMING  PHASE    VARIANCE 


COEFS 


PROVIDE 

ZERNIKECOFFFICIENTS 

FOR   POLYNOMIAL    FIT 

Of    MEAN.    TILT.    ANO 

FOCUS 


PHVAR 


BEGIN   ANGLE    ITERATION 
LOOP   FOR   TARGET    PATH 

INITIALIZE    ACCUMULATED 
FLUENCE 

FLUX-0 


figure  3„ 9    Bain  Program  Flow  Diagram  (coct) 
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figure  3»10    Haiu  Program  Flow  Diagram  (coct) 
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Pigur€  3«11    Bain  Program  Flow  Diagram  (coct) . 
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Figure  3.12    Bain  Program  Plow  Diagram  (cont) 
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figure  3„13    Bain  Program  Flow  Diagram  (coct) . 
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figure  3„14    Bain  Program  Plow  Diagram  (coct) 
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figure  3.15        flaia  Program   Flow   Diagram    (cont) 
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figure  3. 16 


fiain  Program   Flow   Diagram    (cont) . 
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33.   SCEiCCllNE  DESCEIfTICN 


1-   ilSCfiB  and  5CM 

as  noted  in  tie  preceding  chapter,  absorp  tron  acd 
scattering  effects  are  treated  indent icaliy.  The  AfSCSE  ard 
£CAT  algorithms  perfciai  the  integration 

r  atm  (q  5i 

-/  a(h)dh  {'J"j) 


T   =   exp 
a       ^ 


h 


T   =   exp 
s 


h 


atm 

a(h)dh 


(3-6) 


Ihe  result  is  the  tr  aisaission  due  to  the  total  a rscrt tioc 
cr  tctal  scattering.  Ihe  correction  for  slant  range  rs 
applied  iithin  the  angle  interval  loop  or  the  main  crograa 
and  is  siatly 


(sec0) 


=   T 


(3.7) 


is  the  zerith  angle  of  the  target  with  respect  tc  the 
transnitter . 

Ihe  extinction  coefficients  for  absorption  (a)  arc 
for  scattering  (o)  are  provided  to  the  subroutines  by  a  call 
to  routines  ALFA  and  4IFS  respectively. 
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TABU  I 

AESCEB 

and  SCAT 

Prograu  Variable  Definitions 

Variable 

Fortran 
3uE 

na  He 

~ 3  a  in 

Definition 

1     (MSChB) 
a 

I 

TABSC 

Total  molecular 
transmission 

1S(SCA3) 

I 

TSCA1C 

Total  scattered 
transmission 

h   (Eeters) 

atm 

HATMC 

HATMC 

Heignt  of 
atmosphere 

h  .  (ff€t€IS) 

HT 

HTfiAtiS 

Height  above  MSI 
of  transmitter 

a(hj  (Jin)"1 

ALP 

— 

Absorption  coef.       , 
at  given  altitude 

a<t)  (kn)  1 

A1S, 

- 

Scattering  coef 
at  given  altitude 

- 

2J 

N2 

#  of  integration 
intervals 

2-  III-   and  A1FJ 

Ihese  twc  routines  are  also  identical.  They  provide 
the  extinction  coefficients  tc  ABSOEB  and  SCAT  for  a  speci- 
fied altitude.  At  each  altitude  where  a  coefficient  value  is 
desired,  a  linear  interpolation  is  performed  between  data 
points  supplied  hy  tie  user.  Data  statements  precede  each  of 
these  subpr cgrams,  and  it  is  with  these  statements  that  the 
absorption  and  scattering  data  should  be  entered.  Units  for 
the  coefficient  and  corresponding  altitude  should  be  km   and 
km  respectively. 


63 


I 

TABU 

II 

AXIS  and 

A1FA  frog. 

cam  Variabl 

es  and  Definitions 

Variatie 

Fortran 
III 

name 

r-     ■ 

Definition 

0(h)  l**)"1 

S 

AIS 

Scattering  ~oef. 
at  specified 
altitude 

att)     (ka)"1 

A 

Aif 

Absorption  coef. 
at  specified 
altitude 

h  (meters) 

H 

- 

specific  altitude 

All  (M) 

— 

Altitude  data 
list 

— 

ATA  (Kl) 

- 

Absorption  data 
list 

ATS  <U) 

— 

Scattering  data 
list 

NL 

i-Jumter  of  points 

in  data  list          4 

I 
i 
J 

3 .   C OC £1 

Jcr  a  gaussiar  team  with  constant  phase,  the  ititial 
amplitude  distribution  at  the  transmitter  aperture  is 


U(r)  =  U0exp 


-(r/w] 


(3.8) 


where  UQis  the  amplitude  and  w  is  the  spot  size  JEef.  21]. 
Ihe  purpose  of  UOCST  is  to  compute  the  constaot  U0  f  or  a 
given  aperture  power. 

Ey  the  scalar  wave  approximation,  the  intersity 
distribution  is  [Bef.  22] 
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I(r)  =  |U(r 


(3-9) 


3o  relatt  tie  field  distribution  to  the  power,  toe  iriter.sity 
distr ilu tier  is  integrated  over  tne  aperture. 


P,  = 


*^r . 

1 


I (r) dr 


(3.  10) 


£ubs titu tire  3.8  and  3.9  into  2.10 


P.  = 


i 


U0  exp 


-(r/w) 


dr 


{3.  11) 


then  rearranging  teras,  produces  the  expression  f  or  U 


U   = 


A*f 


p.       iexp 

r 


-(r/w) 


2  dr2 


(3.  12) 


2he    intecration    limits    r      and    r.    are    the   radius   of    the 

o  1 

transnitter    and    the    crscura tier,    respectively. 

10C£T   perioms    the    integration   in    3.12    using    the 
trap€2oical   rule.       A   call  is    made    to    subroutine   FIE1L    tc 
evaluate   exp  [-  (r/w)  2]    which    is    siiaply    U  (r)     with    U0=    1.    for 
this    reason,   U  0is   defined   as    ucity   when    FIELD  is    called   r_y 
U0CS1. 
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TAEIE  III 
DOCST  Erogiaa  Variables  and  Definitions 


V  5riaxl€ 

Fortran 

3ul 

na  He 

~Hain 

I  o    (*)  2* 

UO 

- 

ft     (Katts) 

P 

P1C1AX 

- 

N 

ill 

ro(ff€t€IS) 

HO 

- 

r  j_  (mete  rs) 

EI 

- 

^2     " 

D(i)     (  w  }    ni 

i 
US 

- 

*     (mete  is) 

BMBAI 

- 

Definition 

Normalization 

constant 

Total  aperture 
exit  power 

#  of  Integra tioc 

increments 

Outsr  radius  cf 
transmitter 

Radius  of 
onscurat ion 

Field  ampiit  ude 
at  radius  r 

Radius  of  sp  ot 
at  aperture 


1  *  11112 

IlilD  calculates  tie  field  distribution  for  a  axi 
symmetric  Gaussian  beam  with  constant  phase. 


U(r)  =  U0exp 


-(r/w) 


(3-  13) 
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TABU    IV 
FIELD    Prcgraa    Variable   Definitions 


Variable 


-i 


Fortran   name 
SuB  "lain 


D  dr)  <fc)  m  l     UE 


<fc)  2c  -i   UO 


w  (meters)    3MF.AI 


Definition 

Field  amplit ude 
at  r 

Normalizatio  n 
constant 

Fadius  of  sp ot 
at  aperture 


5  .   iJE FID 

In    the    Fraunbcfer   regicn,    the    amplitude   1  istritutici. 
can    be    fcund  by    takinc    the    Fourier   transform   of   the    aperture 
cistr ituticr  [Ref.    23].      Using    pclar    coordinates    and   noting 
axial    s^iimetry,     0*(p)    in    the    far-field    can    be    expressed    as 


U(p)  = 


Az 


b       2tt 
J      /~U(r)exp 


0 


i2r 

rz»-rp 


rdrdG 


(3.14) 


Ising    integral    relaticn   for    the    zero    order    Bessel    function, 
eguaticn    3.14   can   be    simplified    to   £Bef.    24] 


U(p)     = 


2tt      f 
Az    J 


U(r)    J 


2iTpr 
Az 


rdr 


(3.  15) 


67 


Ihe    irter.sity   in    the    iar-fiela    is    then    given    ry   I  (r)  =  |U(r)l 
cr    in    terns    of    eguaticn    3-15 


Kp)    = 


2tt 
Az 


)T/' 


U(r)  J 


2irpr 


rdr 


(3.  16} 


where    i    and   a    are    the   radius    ci    the    transmitter   ind    the 
radius    ci    the   obscuration,    respectively.      [Eef-    25] 

I^BIID    evaluates   3.16    at    a   specified    number    ci 
increments    in    the   iar-rield    and    assigns    these   values    tc   an 
arra}    I  (i) .    F(i)    is    rcrnalized    with    the    total   aperture    fower 
so    that    the    HTF    produced   would    he    unity    at    the   origir.     Kcte, 
however,    that   when    tie    MTF    is    commuted    in    this   case,    the    IMF 
is   normalized   so  as    tc   produce    an   Irei    {intensity    relative} 
value    wren    integrated. 


M     Cp) 
a 


2j.z  : 

(Az)     I 


'/ 


Kp)J 


27Tpr\       — 


Az 


pdp 


(3.  17) 


I0  is    deterained  by    evaluating    eguation    3.16    with    r=C.O- 
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TABLE    V 

. 

I&BFLd 

Program 

Variables 

and   Definitions 

Variatle 

Fortran 
111 

nane 
~Haj.ii 

Definition 

k(rn)"1 

CKR 

- 

(2ir/A) 

P 

E1 

- 

fiadius    in 
far-field 

t 

EO 

- 

Radius   in 
aperture 

I0  (Katts) 

PMA2C 

PUAXG 

Far-field  0Q-axis 
intensity    for    a 
diffraction 
limited    ream 

1(L)       (ft/flj2) 

F(I) 

TISC 

Irradiance    in 
far-field 

Ha(p) 

Gil) 

IEEMTF 

Aperture    MTF 

DX2 

DX 

Aperture   MTF 
increments 

Et(Ratts) 

P 

P1C1AL 

Total    power    in 
the   aperture 

€  -  I£i£EL 

If  the  'times  diffraction  limited  number*  (TDF1M1)  is 
not  specified,  this  routine  computes  the  Irel  value  due  tc 
team  gualit^y  which  is  then  used  to  compute  TDFLMT  in  the 
main  program.  The  calculation  performed  is 

Irel  =  2tt  /  M^  (~p)  Mfo  ("pj  ^d^  (3.18) 

where  Mais  the  aperture  MTF  and  Mils  the  MTF  for  the  team 
guality  phase  screen.   Ma  is  defined  in  subroutine  FABFID  so 
that  this  integration  produces  an  Irel  value.  Mj-,  is  the  team 
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guality    lf.ll    as    computed    by    surrcutrne    M1F3Q.    1DFL  Ml    is 
commuted   in    the    mair.    program   arid   is    1/   Irel    . 


TABU    VI 

EQXREL  Prociai  Variatles  and   Definitions 

llL±3±ll  F2r^xan   name 

3uZ            HTaTn  2§ii^iiion 

K                               -                     TEFLtfl  limes    dixf. 

limited    # 

Il€l  REL  TE£ 

M     ( p)                     All)             IEEM1E  Aperture    MTF 
a 

Mv.(Pj                      F{I)              -  MTF    for    phase 

D  screen 

N4                  N4  #    of    iteraticrs 

for    M1F    calculaticn. 


7 .   EjljAV 

If  the  EMS  phase  distortion  parameter  (WAVEEQ  is 
rot  specified,  this  xcutine  computes  it  based  on  the 
approximation 

Irel      =      exp    (-a2)  (3.  19) 

2tt6 

*here     a   =   ? -and     6  is    the    EMS    value    cf    the    phase 

A         rms  r 

distcrticn  at  the  aperture.  Irel  as  the  intensity  degrada- 
tion due  to  neam  quality  and  is  egual  to  1/(N)2-  N  is  the 
input  parameter  IDFLMl.   DE1MAV  first  evaluates  a2     using 
eguaticn  3.19  and  ther  uses  this  value  as  a  starting  pcint 
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to   centum    a   acre    accurate 


by   iteration    using    subroutine 


EQIRII. 


is  used  tc  commute  the  beam  quality  M1F  and  tie 


learn  cuality  Irel  as  in  equation  3.18  .  a2  is  adjusted,  and 
tie  process  repeated  until  the  Irel  found  by  equation  3.1S 
is  egual  tc  that  determined  by  1/(N)2. 


TABLE  VII 

EFTHAV 

Program 

Variables 

and  Definitions 

V 

ariatlt 

"2tt6 

rrnc 

2 

Fortran 

511 

VASEw 

Name 

"lain 

7  A  5  B  C 

Definition 
Constant 

X 

6 

rms 

WAVEEQ 

WAVF.EQ 

Phase 
distortion 

A 

parameter 

1/N2 

EE10 

V(IDflMT)  2 

6.      ffIFE£    and    BM 

Ihese    tvc   routines    calculate    the    beam    quality    M1F 
array-    lie    tlF    is  [Bef.    26] 


M,   (p)    =   exp 

D 


(-k2tz  -  cM) 


[3.20) 


where  Ch<P)  is  the  autocorrelation  function  of  the  phase  and 
is  tie  pbas€  variance.  CJp)  is  assumed  to 'be  Gaussian. 
letting  C^p)  •=  azexp[-  ( p/L)  2]  f    where  L  is  the  phase 
correlation  length,  results  in  [Ref.  27] 


71 


M  (p)  =  exp 


2tt6 


rms 


1-exp 


-(o/L) 


(3.2}) 


rms 
— t—  is  the  wavelength  HAS  chase  distortion  and  is  ar 

input  parameter  (WAVIEQ)  .   I  ,3150  a  user  input  (SCAIBC)* 

will  default  to  1/5  tie  diameter  of  the  aperture  if  ret 

otherwise  specified. 


1ABLE  VIII 
M1FEQ  and  BM  iiogram  Variables  and  Definitions 


VarijlJ 

~2u6 

rms 

j 

2 

Fortran 
III 

VA3EC 

F(I) 

N4 

nace 
~ Hain 

VAEEg 

BCfllf 
N4 

Definition 
Beau    quality 

A 

ebH3 

variance 

Beam    quality 

air 

#    of    KTF 

k  (m)   ' 


integration 
increments 

(  2  7T/A  ) 
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S.      ISCXEB 

ISC1EB    calculates   the    zenith,   iso^ianatic    angle    for 
use    in    the    adaptive    cities    portion   of    the    prograa .    The    ancle 
is    given   by   £Bef.    28] 


Go    =    .314 


2.91 
6.88 


sec 


atm 

/"c'(h) 


(h)     h    3dh 


(3-22) 


Por    th€    ceai   zenith    case,     sec6=l   .    c2    is    the    refractive 


n 


index  structure  constant  and  is  calculated  hy  subroutine 
CN2H  as  function  of  altitude.  See  Fried  £Bef.  22]  fcr  a 
discussicn    cf    the   ancve    angle- 


TAB1I    IX 
ISOTfiB  Program   Variables   and   Definitions 


Variables 
Oo    (rad) 


/, 


C2     (fl)      ' 

n 


9    (dec) 


ouB  lam 


ISOAKG 


CN2 


ISCASO 


N3 


OMEGA 


Definition 

Zenitn 

isoplanatic 

angle 

Befractive  index 

structure 

constant 

Integration 
intervals  for 
turbulence 

Zenith  angle 


1C-  JCIKT  and  SIK1IS 

If  full  adaptive  optics  are  used,  amplitude  scintil- 
laticn  eifects  will  act  to  degrade  the  performance  of  the 
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con venticnal   phase    ccirecting    adaptive    oj.ti.cs   system.  This 

subrcutire    computes    the    log    amplitude    variance   of    the  scin- 

tillaticts    xor    zenith    trans aissicn-    The    approximation  used 
is 


g2    =.56k 
z 


7    yh 

'/$f    a  tip  5/ 

h    C2 (h)dh 
n 


(3.23) 


For  the  eff  zenith  case 


2  =  a2  sec(9) 


(3. 24; 


vhere    0   is    the    zenith    angle    [Bex.    30  J. 
ance    less    due    to   sciitiliat ion    is 


in. 


relat i  ve  ir  raci- 


AMP  LOSS 


exp( -o_sec  (9 


5) 


The  correction  fcr  tte  off  zenith  case  is  applied  ii  sulrcu- 
tine  2INI1S  in  the  ancle  loop  cr"  the  main  program  when  full 
adaptive  cities  are  specified.  This  loss  factor  has  teen 
shown  tc  he  limited  tc  approximately  .55,  therefore,  a 
default  value  is  used  in  the  main  program  in  tne    case  where 
eguaticr  3-25  yields  a  value  less  than  0.5  £Ref.  31]. 
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SCIMIS   and    SCINl   Program   variables   and   Definitions 

^ilij^ii  Fortran   Nane 

SuE  '  HaxJi  Definition 

k     (a)"1  CK  -  2tt/A 

C2     (ir)    3  CN2E  -  fiefractive  index 

|  n  structure 


constant 


a2  SIG52           SIGXZ                Log    amplitude 

1x  variation 

iUp  less             TAtll             TAMP                   loss    caused    tv 

]  scintillation 


i 


«. 


11.     CN2£ 


CN2E   computes   the   index    of    refraction   structure 
constant   as   a   function   of  altitude.    C2    is    determined    usinc 
Hufnacei's    iiodel    with    an   added    term   to   include    the    near 
surface    effects.      £  fief  .    32  J    £Bef.    33] 

C2    =     (5.94xlO--)h-exp(T^)  (3.  26) 


+    (2.7xl0_1  6)exp' 


1500 


♦    (lO-'^exp^) 


Here,  the  units  of  h  are  meters. 

12.  JHCTEE 

Ihis  subroutine    computes    Yura's   [fief.    34]]    lateral 
coherence    length#Po    ,    in    terms    of   Fned's    [fief.    35]    coher- 
ence   cianeter,    r». 
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Po    =    r 0/2 . 1 


(3.21) 


Substituting   Fried's    definition    for   r 


Po    = 


2.1 


?    Ql  f    atm 

|^|i-  k2  /      C    (h)dh 
6.88  /,        n 


_3 


V* 


(5.26) 


TABU    XI 
EHOTRB   Procram   Variables   and   Definitions 


Variable 


k    (m) 


-  1 


-2 


C2      (B) 

n 


-/■> 


Po   (E€teiS) 


Fortiaii    11  a  me 
III  "I^in 

CK 

CN2 


BHC 


iJECO 
N3 


2tt/  A 

Refractive  index 

structure 

constant 

Coherence   length 

#    of    integration 
intervals   ror 
turbulence 


13-    *jll   and    ETFJIl 

JI3    AND    MTFJ1I   function    aigorithmically  the    san€    as 
Eft    and    MIIIQ.      Tcgetler,    thev    confute    the    jitter    HTF    anaj. 
HIFJII    is    called  by    JIT   for    each    radial    increment.    The 
jitter    Mil    is   given   tv 


2tt^ 


M,  ( p)    =   exp/k^(2V 

-J  \  O 


(3.29) 
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where    2q     is  a    user    specified    input.    As    with    fceaa     guaiity, 
this    cuartity   is   applied   ac    the    aperture. 


TA£IZ    XII 
C11   and   aiFJIT    Urogram    Variables   and    Definii.ticES 


Variable 


k    (a) 


r  1 


Fortran    Name 
Sub"  lam 


CK 

F 


JIIM1F 


^iiiiiticn 

2tt/A 
Jitter    MIF 


1^4-     IB1EHS 


Ihis  subroutine   calculates    tiiat   part    of   thernal 

tloctting    net   dependent    on    the    slant   patn    of    the   beam.    Ihe 

eguaticn    used   by   GUISAVG  for   the    phase   distortion    due    tc 

tnernal    ilccming   is 

2jt    n-1    y-1    fX'    I(x',y)dx'l 
A       Po       Y     Jm  J 


A(t)(x,y)  = 


X 


atm 


h 


a (h) exp  (- 


-sec8/ 


a(h)+a(h)    dh 


dh 


(3.20) 


't  V0cos^    +    toh 

Note    that    tie   first    term  in    this    expression   is  invariant 
fcith    respect  to    team    rath   while   the   second    term   contains 
such    catt    dependent    variables    as    wind    ,    wind    due    tc    slew, 
and    eitiiction    coefficients.    The    first    term   is   evaluated    fy 
this    sutxoutine    and    the   remainder    evaluated    later    in    the 
progran    inside    the    sheet  angle    iteration   loop   by    subroutine 
AV. 
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.Standard  values  used  icr  n  -1,  p0/  and  "Y  are 
2.72x1C4,  1.01x1Cs  J/a3/  and  1.4,  respectively,  k     phase 
screen,  i(i,j)  ,    is  cc:s true  ted  by  iteration  ox  the    irtersity 
integral.  Ncte  that  picgram  performs  tnis  integration  even  a 
naif  place  cf  the  aperture  defined  by  -t<y<fc  and  0<x<r  where 
t  is  the  aperture  dianeter.  Ihe  results  are  adjusted  later 
in  s  cir  "entire  PHVAR  fcx  this  metned. 


TAE1E 

XIII 

PFTPHS 

Program 

Variables 

and  Definitioxs 

Variable 

Fortran 
III 

Na  me 
~ lain 

Definition 

k  (u)"1 

CK 

- 

2  7T/A 

».-' 

const. 

*~ 

Refractive 

index  term 
(2.  72x10*) 

Y 

const. 

- 

Ratio  of  spa  cif ic 
heats  (1.4) 

Po 

const. 

— 

Atmos.  pressure 
(1.01x10s)  N/m2 

*  U/Y) 

PH(i,j) 

PH(i 

rj) 

Inermal  hlooming 
phase  array 

15.  BESFHS  and  CCLFFS 

Cnce  the  phase  screen  has  been  constructed  by 
subroutine  fETPHS,  subroutine  BESFHS  removes  the  near,  tilt, 
and  fects  curvature.  Ihe   result  is  tne  phase  aberration  due 
to  theriral  blooming  alone.  Zernike  polynomials  are  utilized 
for  expressing  these  phase  distortions  with  subroutine 
COEFFS  providing  the  required  ccexf icients-  These  pclync- 
nials  are 
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Z0  (x,y)     =   a0                                     {mean)  (2.  21) 

2i   (x,y)     =    ai  x                                (x    tilt)  (3.22) 

Z2  (x,j)     =   a2  y                              (y    tilt)  [3.23) 

23   (x,y)     =   a3   (x2+y2)  +  at      (focus)  (3.j4) 

Ihe    expaisicn   coefficients    axe    given    fcelow    and   are    ccmp uted 
xelative    to   a    uniform    aperture    weighting    function    *(x#y). 
Mx,y)    =    1    inside    the  aperture   and   zero   elsewheie.    Ihe    inte- 
grals   have    reen    nultiplied    ty    2    to    compensate    for     the    half 
plane    integration   of    subroutine    PHHPHS. 

f0=    2a0       $(x,y)W (x,y) dxdy  (3.25) 

fi=    2ai       (j)  (x,y)  W  (x,y)  xdxdy  (3.26) 

f2=    2a2       <|>(x,2)Ji  (x,y)  ydidy  (3.37) 

f   =    2         (a    (x--»y2)+a    )     (x,y)  W(x,y)dxdy  (3.28) 

3  3  4 

Ihe  final  chase  correction  is  given  by 

<Kx,y)=<Mx,y)  -f  iZ  i  (x,y)  -f  2Z2  (x,y)  -f  aZ3  (x,y) 

-f 0Z  0 (x,y) 


(3.3S) 


lor    a    discussion   of    least  sguares    fitting    of    Zemike    pclyrc- 
irials,    see    £Ref.    36]. 
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16-    iB\£* 

Tre   phase   variance    due   tc    thermal   blocmia  g    is 
commuted    iy    this   subicutine    using    the    residual   phase    screen 
provided    ry    RESPES.     Kcraaiized    with    respect    to    the    apeituic 
field,    tie    variance    is    given   iy 


oz    =    U(x,y)  (j)2(x,y)     +      U(x,y)  4)(x,y)  (3.40> 

U  ( x  ,  y ) 


TA£IE    XIV 

j 

i 

\ 

PHVAB,    SESPHS 

#    and  CCZ.FS    Program 

Variable    Definitions 

j 

i 

Variable 

Fortian 

5U — 

Name 
~Hain 

Definition 

<J>(x,;y) 

PH{i,jj 

PH  (l/j) 

Pnase    array 

i 

SIGSQO 

S1GSC.0 

Pnase    variaace 

1 

ao"    a. 

Al-    A4 

- 

Zernike    coeffs. 

Vi> 

PHMIAN 

phtitx 

PHTIll" 
PHI CCUS 

- 

The    expansion 
coefficients 
relative    to    a 
uniform    weighting 
function 

1 
i 

- 

J 
_J 

17.    iV 

AV    evaluates  the  path  dependent  part  (second  tern,) 
cf  the  tiernai  blooming  equation  3.30  .  (see  section  PRlPhS) 
a(h)  and  ®  (I)    are  the  absorption  coefficient  and  scattering 
coefficient,  respectively.  V0  is  the  user  specified  wicd  and  co 
is  the  slew  rate.  The  denominator  represents  the  total  tran- 
sverse i»ind  component  across  the  'team.  0  is  the  zenith  angle 
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and  5  is  the  angle  cf  attack  cf  V0  with  respect  to  the  teak. 
V  has  teen  assumed  to  re  opposite  in  direction  and  parallel 
to  the  target  to  trarsmitier  1CS  slew  motion.   Note  alsc 
tiaat  V   sill  be  applied  as  a  ccistant  the  entire  length  o.f 
the  team  path.  This  it  discussed  in  chapter  one. 


TABU    XV 

AV  Prograa  Variables  and  Definitions  J 

Variance     ^ortzdn   Name  | 
5uE   ~   Ills        definition 

a{h)  +  a(h)    ALP      -          Total  extinction  \ 

coefficient 

a\h)     (J<m)    A1S      -          Absorption  coeif.  j 

sec  (9)        SECCV.G        SECCKG  secant  or  zenith 

angle  j 

l0(m/sec)     VO       VO         Kind 

PHIICT    PHILCT      Slew  rate  j 

I 

h          (meters)     HAT2C           HA1MC                 xieight    of    ttfe  j 

atm                                                                       atmosphere  j 

h       {aeters)       HT                  H1RAK5             Height    of   the  i 

t                                                                            transmitter  j 


16.  E1CCM 

Ihe  purpose  cf  BLOOM  is  to  provide  reasonable 
results  for  thernal  tlccming  degradation  when  the  Irel 
values  aie  below  approximately  0.3.  In  this  region,  the 
exponential  Strehl  relation  predicts  unacceptabiiy  severe 
results,  therefore,  *len  the  phase  variance  is  greater  than 
1.2,  ELCCM  computes  ar  Irel  value  based  on  GUISMTF  results. 
GUTSETF  is  a  full  wave  optics  propagation  code  using  FFTs. 
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Curve    nt    ccly  nonials    were    developed    using    GUlSt-lIF    results 
for    lets   a    uniform    arc   a   truncated   Gaussian   aperture    distri- 
iuticn.      !he  resulting   pciy Hernials   are 

irel    =      1 (3.41) 

-.08705    +    2.9148a        +    .1723a2 


Irel    = 


1.2877    -    2.6491a    +    4.09603a2  (--42) 


If    the    profile    under    consideration    is    a    truncate!    Gaussiar, 
i.e.    the    aperture    diaceter    greater    than    trie    waist    diameter 
ex    the    ream,   then    the   routine    uses    the    truncated    Gaussian 
polynciiai.    Otherwise/    a   comhination    of    the    two    are    used, 

lreltb    =       j~l-(d/b)  21  lrelu    +    (d/b)2Irelt  (3-WJ 

kuere    d    is    the    apertue    dianeter    and   b    is    the    waist    cianeter 
cf   th€    teas.      [Befo    27] 
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TAEIZ  XVI 
EIOOM  Program  Variables  and  Definitions 


Variable 

a2 
II6ltb 

II€lu 

IX€ltg 


Fort i an  Name 
Sub"     ~~ Ham 


S 

T 

TRU 
TfiG 


SIGSQ 

TEICCM 


Defin 

Phase 

Iotal 
prcdu 

Curve 
for  u 
apert 

Curve 
for  G 

apert 


ition 

variaa  ce 

bloomi  ng 
ced  Irel 

fit  Irel 
niform 
ure  dist. 

fit  Irel 
aussian 
ure  dist. 


1S.    gllJIM 

EIJATM    computes    the    atnospheric    HTF    and   if    speci- 
fied,   a^flies   a    tilt    due   to    turbulence   correction.    Ihe    Elf 
is    giver   by 


Mt(p)    =    exp 


-(p/D)  Vj    [l-d(p/D)  1/3J  (D/p0) 


v, 


(3.  iii*) 


fchere    C    is    the    ajertire    diameter,  p0  is    the   coherence    diam- 
eter   and   d    is   defined   as    1-A£A£.    ADAP    is   the    fractional 
residual   tilt   due    to    turbulence   and   is   a    user   specified 
parameter.    If   ADAP=0,   M      will   represent    the   fully    turbulence 
tilt    corrected    MIF.    If   AD/P=1,    M      will   be    the   uncorrected 
MTF    aid    uav    be    writter    as 


Mfc(p)    =   exp|-(p/p0) 


v 


(3.45) 


ALAP  nay  te  any  value  between  0  and  1 


TA£IE    XVII 

H1FATM 

Program 

Variables 

and   Def initioms 

Variable 

Fortran 
5u5 

Ha  me 
~flam 

Def rait ion 

Mt("p) 

F 

F3 

Atmospheric    'All 

p0  (a  etc  is) 

RHO 

HBO 

Coherence   dl  a. 

E     (meters) 

D 

DIA 

Aperture    dia . 

i 

2C-    IBEFEL 

I£E£EL    computes    the    a^frcxiMce    Irel    for    a    tilt 
corrected    system.       This    value    is    used    for    comparison    with 
the    Irel   produced  by   a   prefect   adaptive    optics   compensated 
systeir.    Ihe    purpcse    of   this    ccirparison   is    to    determine    if 
noise    will    degrade    tie    AO   system    to    an    extent    as    to    sake    a 
full    AC    system    undesirable.    If    thrs   is    the    case,    the    program 
sill    apply    tilt    correction    only    to    the    team. 

21-    OD£flO 

Subroutine  FKERHO  calculates  the  coherence  dianeter 
that  *culd  result  in  the  Irel  value  achieved  by  a  perfect  AO 
systen.   Ihis  .is  done  by  successive  calls  to  subroutine 


IRBBEI  and  iterating  p 


This  new  Pn  then  becomes  a  factor 


u  o 

in  the  AC  compensation.  Specifically  it  will  re  used  by 

subrcttine  flTFATM  to  produce  the  atmospheric  MTF  of  the  AC 
corrected  system. 
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TABIE  IVIII 
1EBBII  and  FKDBHC  Program  Variables  and  Definitions 


Variabl 

€ 

For tian 
III 

Name 
ETaln 

Ka(P) 

F2 
F1 

IBBMIF 

Iiel 

EXD 

EXD 

Atmospheric  M2F 
Aperture  SIP 


Irel  or  perfect 
AO  system 


2  2-  HI  ISO 

She  call  to  HUSO  is  nade  when  adaptive  optics  have 
rot  teen  user  specified  and  iscplanatic  calculations  have 
Let  been  inhibited.  Also, if  the  signal  to  noise  ratic  at  the 
AO  sensor  is  such  that  the  AO  system  will  provide  tilt 
correction  cnly,  I1TIBC  is  called  ,  again,  provided  iscpla- 
natic  calcuiations  have  not  been  inhibited. 

Ihe   purpose  cf  T1TISC  is  to  include  beam  wander  due 
to  iscplaEa tism.   The  2-sigma-p  tilt  is  computed  and 
combined  with  the  jitter  2-sigma-p.   This  comtined  tern  is 
then  used  tc  compute  the  jitter  MIF. 

23.  5EIIOT 

EEIIGT    simply    averages    the   relative    intensities    due 
to    a    aultiplica tive    aid    an    RSS  (root    sum    squared)     approach    to 
combining    thermal   blccming    and    the    other   propagation 
effects. 


T     1       Irel        +   Irelm 

Irel  =      rss        m 


(3.46) 


Irel   is  the  result  ci  the  multiplicative  approach 
m 


Irel   =  Irel.,  *  Irel 
m        tb        o 


(3.141) 


where  Irel,bis  the  Ii€l  due  tc  thermal  hlocming  and  IrelQ  is 
the  Irel  due  tc  team  cuality,  jitter,,  turbulence  ,  isc^Iar.a- 
tism    and    adaptive    optics   effects-    I  he    HSS    abroach    is    giver. 


Irel 


rss 


1+ 


1 


Irel 


"    -    1      + 


1 


Irel 


-    1 


tb 


(5.48) 


TAB1E    XIX 
EilTOT   Program   ¥ariafles   and  Definitions 


Variable 

Fortran 
III 

Na  me 
~Harn 

Definition 

Iielt 

T 

T101M 

Total  Irel 

Il€lrss 

TR 

TBSS 

Irel  by  ESS 
method 

irel 


m 


Irel 


Irel 


tb 


TM 


EXD 


TMUI3  Irel   by 

multipljLcatl  ve 
method 

EXD  Irel    of    perfect 

AG    system 

T1RUE  Irel    due   to 

a-i.1    other   effects 

TEICCM  Irel    due    to 

thermal    blooming 
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24.     ISOPLA 

This  routine    calculates    the    MTF    that    characterizes 
the    isoplanatic    effect    on    the    predictive    or    "look-ahead" 
adaptive    optics    system.       The   look-ahead    angle    is    the    major 
input    parameter    to    this    routine.       The    Fortran    code    was 
written   by    D-L.     Fried  [Bef.     38].       Fried    developes    the 
isoplanatic   dependency   of   the    AO    system    in    terms    of    the 
effective   antenna   gain  of   the   laser    transmitter.    The   MTF 
that    Fried    formulates   is  given   by 


M.         Cp) 
ISO 


■/ 


-2  -5 


144.88-X 


-*v> 


■  sec  (9)  ■ u(vsec(6/p)  , 


(3.49) 


where 


x    =    vsec ( 0/p) 


(3.50) 


and 


U(x,<j>)     =    /C2<(l+(xh)        -h 


1+2  (xh)coscf>+    (xh) 


-h  [l-2(xh 


) cos    +     (xh)  l 


■dh 


(3.51) 


0   is  the  zenith  angle  and  v  is  the  target  lead  angle. 

For  detailed  discussion  of  the  theory  and  explicit  develop- 
ment of  the  Fortran  code,  see  [Ref.  39]. 
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JO    commutes    tie    zero    order    Bessel    function    based    en 
the    infut   argument.       Ihis   routine   is    called    ny   FAEFLI    it    the 
calculation    of    the    fai-fisld   irradiance. 
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APPENDIX  A 


GUTSAVG  INPUT  FILE 


INPDl    DATA     EI1I; 


LASEE: 

CLIMATE: 

CASE: 

DIA 

DIAOES 

3EAMSZ 

HAVE 

PTOTAI 

THSEE 

HGBNE 

TDFLKT 

UAVEEC 

SCALEC 

HTEANS 

HSAT 

THETKX 

LOFF 

SHOO 

VO 

SIGJI1 

ADAP 

AOBICfl 

AVGSPT 

NFLAGC 

NFLAGA 

NOISC 

ABSLCZ 

XJT 

BilLTE 

NA 

N1 

N2 

H3 

H4 

H5 


CO    CH    EDL  5P(9}     TBANSITIGN 

BIL-LATITDDE    SUMMEE,     CLEAH    DAI 

1 


TILESCCEE    DIAMETEE 

CEKTEAL    OBSCCBAT1CN     DIAMETEB 

GACSSIAh    WAIST    TfcEU    I2LESCCPE 

CSVITY     WAVELENGTH 

AEEBTDBE    TOTAL    PC'nEB 

TCEEDIENCE    SEEING 

EE1GEI     CF    GBCUND    AECVE     MSL 

TIMES    DIFFEACTION    LIMITED 

E£S     VAVIS    DISTCETICN 

PEASE    CCBEELATION    LENGTH 

TEAfcSMlTTEB     fcEIGEl    ABOVE    HSL 

SAIEILIIE    ALTITDCE 

HAXIEUJ    ZENITH     AfcGLE 

IIIGE1    PATH    CFFSET 

IDEA'S    1UBB.     COHERENCE     CIA 

SIND    VELOCITY 

2-SIGMA-P    JITTEB 

lllC.    EESID-    TUBE.    TUT 

E1CCMING    COEBECTICN     (0-1) 

SECT    AVEBAGINF     FACTCF.     (0-1) 

IKCLLDE    CLOUDS     MCCEI 

OSE     FDLL    ZONAL    A-C    SYSTEM 

INKIEIT     ISOPLANA1IC    CALCS 

ZENITH    IRANS    FOB     i-C    SENSCE 

TJBGET    BAOIANT    INTENSITY 

A-C    SYSTEM    BANDKIETH 

NDHEEB    CF    A-C    SYSTEM    ACTDATOBS 

KCCEIE    CF    ANGLE     l^EEVAIS 

AESCEPTION     INTEGEAIIOK     INTEBVALS 

TIJBEULENCE     INTEGBAI1GN     INTFHVALS 

HTI     1NTEGBAT10N     IKTEBVA1S 

TEEE5A1    BLOCKING    INTEGEATICN    INTEBVALS 


METEES 

METEFS 

METEES 

METEES 

WATTS 

ABC-SEC 

METEES 


METEES 

METEES 

METEBS 

DEGBEES 

METEES 

METEES 

M/SEC 

BADIANS 


W/STEB 
HEBTZ 


4.SS2  1C    tlCECNS 


C.  15C000E+01 
C.  150000E  +  00 
C.  10000GE+02 
0.  M95210E-05 
C. 20C0O0E+O7 
0.  G 

0-300000E+03 
0-  120000E+01 

c  c 

0^300000E+00 
C.3JQ000E+03 
C.  100000E+07 
C.30COOOE+02 
0.  0 

c.c 
c.eocoooi+01 


1 

C.75COOOE+00 
0. 2500001+03 
C.S00OO0E+03 
0-  102400E+04 


0-  102400E+04 
30 
200 
100 
200 
100 
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APPENDIX  B 


GUTSAVG  OUTPUT  FILE 


LASEE:  CC    CB    EDI  5P(9)      IEANSITICN 

CLIMATE;  MID- LA TIT  IDE    SOMMEE,     CIEAB    DAY 

CASE:  1 

DIA  TELESCCPE     EIAMEIEE 

EIAOES  CENIEA1    OESCOEAIICN    DIAMETIf 

EEAMSZ  GAOSSIAN     KAIST    IBEU    TELESCOPE 

HAVE  CAVIIi    BAVE1ENGIE 

PTOTiL  APEBTUEE    TCTAL    FOfcEE 

1HSEE  TUEBUIINCE     SEE1KG 

BGBNE  HEIGHT    OP     GECUflE     ABOVE    MSL 

TDFLHI  TIMES    EIFFEACTICk    LIMITED 

BAVEEC.  RMS     WAVES    EISICEI10N 

SCALEC  PHASE    COEEEIATIC*.    LENGTH 

HTSANS  TRANSMITTEE     IE3GHT    ABOVE     MSI 

HSAI  SATELIITE    AITITOZE 

IHETHX  fUXIMOt    ZZMIf    ANGLE 

LCFf  FLIGHT    PATE    CEESET 

EHOO  YUEA'S    TUEE.    CCfcEEENCE    DIA 

VO  BIND    VELOCITY 

SIGJIT  2-SIGMA-P    JITIIE 

AEAE  FBAC-     EESIE-     IEEE.    TILT 

AOBICH  BLOCMIhG    CCEEICTICN     (0-1) 

AVGSPT  SPOT     AVEBACIKJ     FACTOE      (0-1) 

HFLAGA  USE    PC  II    2CKA1    A-C    SYSTEM 

NOISC  INHIE1T    ISCELAKATIC    CALCS 

ABSLC2  ZENITH    ISAhS     FCE     A-0     SENSCE 

XJT  TARGET    BACIAKT    INTENSITY 

BSIDIB  A-0    SYSTEM    EANLHIETH 

NA  NOMEEE    OF    A-C    SYSTEM     ACTUATCBS 

N1  NUMBEE    OF    ASGLJ.    JMEBVALS 

N2  ABSORETION     INTEGRATION    INTIFVALS 

$3  TOHEOIINCE     INIEGEATION    INTIEVALS 

N4  MTF    INTEGRATION    IKTEEV4LS 

H5  THEBMAI    B1CCBMING    INTZGEATICN    INTERVALS 


METEES 

MEIEBS 

METEBS 

METEES 

BA  TIS 

ABC-SEC 

METEES 


METEES 

METEES 

MEIEBS 

DSGEEES 

METEES 

METEES 

M/SEC 

RADIANS 


B/SIEE 
HEETZ 


U. 99210    MICRONS 


. 150E+01 

. 15CE+00 

.  100E  +  02 

. U99210E-05 

.200E+07 

. OCCE+00 

. 300E+03 

.  120E  +  C  1 

- OOCE+OO 

. 30CE+00 

. 300E+03 

.  100E  +  07 

. 300E+02 

•QOCE+OO 

. OOCE+OO 

. 060E+02 

.5C0E-C5 

.  1C0E  +  01 

.  100E  +  C1 

.  100E+0  1 

0 

1 

. 75CE+C0 

.250E+03 

. 50CE+03 

1024. 

30 

200 

100 

200 

100 
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PATH     ANALYSIS 

BESOLTS: 

STEP 
HC. 

BANGS 
(KM) 

CLE  W 
(MJAD/SEC) 

CF.EGA 
{LEG} 

TIME 
(SEC) 

ATM05.     TRANSMISSION 
SCAT             ABSCEP 

fiHO 
(CM) 

1. 

939-738 

7.1582 

1.08 

2.  63 

0. 99802 

0.53314 

51.  1 

2. 

1CCQ.045 

7.1541 

1.  16 

5.27 

0.99802 

0. 533C2 

51.0 

3. 

10C0.655 

7.1459 

5.24 

7.90 

0.99802 

0.53276 

51.0 

4. 

10C1.579 

7.1337 

4.31 

1C.53 

0.99802 

0.53242 

51.0 

5. 

10C2.805 

7.1175 

5.39 

13.  16 

0.99601 

0. 53195 

50.9 

6. 

10C4.334 

7.0973 

6.4o 

15.80 

0.99801 

0.53135 

50.9 

7- 

1006.  167 

7.073  2 

7.52 

16.43 

0.99800 

0.53064 

50.8 

8. 

10C6.301 

7  .  04  5  3 

6.59 

21.06 

0. S9800 

0. 52962 

50.3 

9. 

10  10.735 

7.0133 

9.65 

23.70 

0.99799 

0. 52888 

50.7 

10. 

1013.465 

6.9786 

1C.70 

26.33 

0.  59799 

0. 52782 

50.6 

11. 

1016.490 

6.94  0  0 

11.74 

28.96 

0.99798 

0.52665 

50.5 

12. 

1019.806 

6.89  3  1 

12.79 

31.59 

0. 99797 

0.52537 

50.4 

13. 

1022.413 

6  .  85  3  0 

13.82 

34.23 

0.  99796 

0. 52398 

50.2 

14. 

1027.305 

6.80  4  8 

14.85 

36.36 

0.99796 

0.52248 

50.1 

15. 

103  1.481 

6.7537 

15.86 

39.49 

0.99795 

0.52088 

50.0 

16. 

1035.935 

6.6999 

16.87 

42.  13 

0.99794 

0.51917 

49.8 

17. 

1040.666 

6.64  3  5 

17.88 

44.76 

0.99792 

0. 51735 

49.6 

18. 

1045.668 

6.584  7 

16.37 

47.39 

0.99791 

0.51544 

49.5 

19. 

105C.939 

6.5237 

15.85 

50.02 

0.99790 

0.51342 

49.3 

20. 

1056.474 

6.4606 

20.83 

52.  66 

0.99789 

0.51  130 

49.1 

21- 

1062.269 

6.3955 

21.79 

55.29 

0.99787 

0. 50909 

48.9 

22. 

1068.319 

6.3288 

22.75 

57.  92 

0.99786 

0.50679 

46.7 

23. 

1074.62  1 

6.2604 

22.69 

60.56 

0.99785 

0.50439 

48.5 

24. 

1081.169 

6.1907 

24.63 

63.  19. 

0.99783 

0.50191 

48.3 

25. 

1067.959 

6.11 96 

25.55 

65.82 

0.99781 

0.49933 

48.  1 

26. 

10S4.988 

6.0475 

26.46 

68.45 

0.99780 

0.45668 

47.9 

27. 

1102. 25C 

5.97  44 

27.36 

71.09 

0.99778 

0.49393 

47.7 

28. 

11C9.740 

5.90  06 

28.25 

73.72 

0.99776 

0.49111 

47.4 

29. 

1117.455 

5.8260 

2S.  13 

76.35 

0.99774 

0.46821 

47.2 

30. 

1125.388 

5.7510 

30.00 

78.99 

0.99772 

0. 48524 

47.0 
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TRAKSEISSICN 

ANALYSIS    OU 

71  tit 

STIf 

TBANEtl 
DUE    TO 
SPREADING 

SS1CN    C0EEFIC1ENTS 
ABE    LOSS          1EESHAL 
HUH    AO              EICCillNG 

RANGE 
SCALE 

MAX    IRRAD 
.     (K./CH2) 

ELUENCE 
(KJ/CM2) 

1. 

0.  11826 

1  . 00  0  0  0 

0.01522 

0. 9999 

0.5779E-04 

C.  1522E-03 

2. 

0.  11821 

1.00  00  c 

C.0152  1 

0.9993 

0.5770E-04 

C. 304 1 E-03 

3. 

0.  11821 

1.00000 

C.0151  8 

0. 9981 

0.5753E-C4 

C.4555E-03 

4. 

0-  11817 

1  . 00  0  0  c 

0.0151 5 

0. 9963 

0.572bE-04 

C.  6063E-03 

c^ 

C.  11810 

1.00  0  00 

C.0151  1 

0. 9938 

C.5692E-04 

C/562E-03 

6. 

0.  11802 

1.00  00C 

C. 01505 

0.9908 

0.5649E-04 

C.  5049E-03 

7. 

0.  11793 

1.000  0  0 

C. 01499 

0. 9e72 

0.5597S-04 

C.  1052E-02 

8. 

0.  11782 

1.000  0  0 

C.0149  1 

0.  3830 

0.5539E-04 

C.  1193E-02 

S. 

0.  11770 

1.00  00C 

C. 01433 

0. 9783 

0.54722-C4 

C.  1342E-02 

10- 

0.  11756 

1 .  00  0  0  c 

C. 01473 

0.9730 

0.5399E-04 

0.  14S4E-02 

IU 

0.  11740 

1  . 00  0  0  0 

C. 01463 

0.967,! 

0.5320E-04 

C.  1624E-02 

12. 

0.  11724 

1  . 00  0  0  0 

C. 01452 

0.96  10 

0.5234E-04 

C.  1762E-02 

13. 

0.  11705 

1  . 00  0  0  0 

C. 01440 

0.9542 

0.  51422-04 

C.  1898E-02 

HI. 

0.  11686 

1.000  0C 

C. 01427 

0.9470 

0.5046E-04 

C.2030E-02 

15. 

0.  11664 

1  . 00  0  0  0 

CO  141  4 

0.9393 

0.4945E-C4 

C.2161E-02 

16. 

0.  11642 

1.00000 

C.0140C 

0.9313 

0.4839E-04 

0.2233E-02 

17. 

0.  11618 

1.00000 

C. 01385 

0. 9228 

0.4730E-04 

C.2413E-02 

18. 

0.  11593 

1.0000  0 

C.0136S 

0.9140 

0.4618E-04 

0.2534E-02 

1S. 

0.  11566 

1  . 00  0  0  0 

C. 01354 

0.9049 

0. 4504E-04 

C.2653E-02 

20. 

0.  11538 

1.000  0  0 

C. 01337 

0.8954 

0.4387E-04 

0.2768E-02 

21. 

0.  11509 

1-0000  0 

C. 01320 

0. 8357 

0.4268E-04 

C.2881E-02 

22. 

0.  11479 

1  . 00  0  0  0 

0.01303 

0.6757 

0.4148E-04 

0.2990E-02 

23. 

0.  11448 

1.0000  0, 

C. 01286 

0.8654 

0.4028E-04 

C.3096E-02 

24. 

0.  11415 

1  .  00  0  0  0 

C. 01268 

0.8550 

0.3907E-04 

C.2199E-02 

25. 

0.  11381 

1.00000 

C. 01250 

0.8443 

0.37862-04 

C.2298E-02 

26. 

0. 11346 

1  . 00  0  0  0 

C. 01232 

0.8335 

0.3665E-04 

C.3395E-02 

27. 

0.  11310 

1 .  00  0  0  0 

C.  01213 

0.8226 

0.3545E-04 

C.  3488E-02 

28. 

0.  11273 

1.00  0 00 

C. 01195 

0.81  15 

0.3427E-04 

0. 3578E-02 

2S. 

0.  11235 

1.00000 

0.01176 

0.8003 

0.3309E-04 

0.3666E-02 

30. 

0.  11196 

1 . 00  0  0  0 

C. 01158 

0.7391 

0.3193E-04 

C.3750E-02 
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CASE    EATA    ANE    CALCULATED    F ACTOBS : 

T3ANSEITTEE     E1AP.ETEB  =  150.00  CM 

CAVII'i    WAVE1ENGTE  =  4.99210P  OM 

APEE1DEE     ICIil    ECHEB  =  2-0000  MH 

TIMES    DIFIEACTICN    LIMIT       =  1.20 
(WIDE    ANG1I     SCATTERING) 

TRANSMIITEB    ALTITUDE  =  300. CO  M 

SATELLITE    ALTITUDE  =  10C0.C0  KM 

MAX    ZENITH    ANGLE  =  30.00  DEG 

FLIGHT    PATE    OFFSET  =  0.0  KM 

BHOO  =  0.5105E+C2  CM 

HIND    VELOCITY       •  =  6. CO  M/SEC 

HEIGHT    OF     GECGNE  =  300. CO  M 

OPT    SEEING    AI    55COA  =  0. 15E+01  ABC-SEC 

BLCOMING    AIAPIIVE    OPTICS    IACTOB       =    0.1000E+01 
IBEAE1ANCE    ABEA     AVERAGING    EACTOB    =     0.1000E+01 

TELESCOPE    DIMENSIONS: 

OUTER    DIAMETEE  =  150.00  CM 

INNEB    DIAMETEE  =  1 5.0C  CM 

GAUSSIAN    HAISI    EIAMETE3  =  1000.00  CM 

EMS    HAVES     EISTCETION  =    0.1031E+00 

PHASE    COEEIIIAIICN    LENGTH    =    0.3000E+02    CM 

FOB    TEE    FELL     20NAL    A-C    MOEEL: 

ZENITH    TEANE.     AT    A-C    SENSOE  =    C.7500E+-00 

TAEGET    EAE1ANT    INTENSIT1  =    0.2500E+03    W/STEB 

ADAPTIVE    CEIICS    BANLHIDTE  =    0.5000E-1-03    HEBTZ 

NOBEEB    CI     ACTOATOES  =  1024. 

THIS    EUN    - 

IS    TEE     EASIC    CCDE     SJITH    CNLY    SOME     HEASOEE    OF    TILT    CORRECTION 
ANE    WITHCOI    AN    ISOPLANAIJC    MODEL 

NUMBEE    OF     1NIERVAX    STEPS  ICE: 

ANGLLLAE    INIEEVAL  =  30 

ABSCBPTICN    INTEGRATION  =  200 

EHC    CALCULATION  =  100 

H1I    CALCULATION  =  200 

IHEBMAL    ELCCMING  =  100 

ZEHITB    LOG    AMPLITUDE    VAEIANCE    =       0.4905E-02 
BELATIYE    IEEAEIANCE    SEDUCTION    =       0.9951E+00 

ZENITH    LOCK-AfcEAE    ANGLE  =    0.5227E+02    DEAD 

ISOPLANATIC    JIITEB     (2-SIGEA-P)     =     0.2004E+07 

2-SIGMA-P     EEAM    JITTER  =  5.00    DEAD 

TDSBULENCE    JITTEB    REJECTION  =     0.1000E+01 

(BESILUAL    =    INITAL    *     ADAP) 


========    EECPAGATION     BESDLTS    ======= 

INTEGRATED    FLUX    CN    TAEGET  =  0.007499    KJ/CM2 

TOTAL    ILLUMINATION    TIME  =  157.97    SEC 
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APPENDIX  C 


GUTSAVG  PROGRAM  LISTING 


c** 
c 
c 
c 

c 
c 

c 
c 
c 
c 
c 
c 

c 


******************************************************* 

GOTSAVG         NAVA1  PCSTGEAEOATI  SCKOOI  VERSION  1.2 

BODIFICATIOK  1.2 

IATES1  CHANGE  EA1I  22  FEE  84 


*  *  4*** ******** 

* 
* 
* 
* 
* 

* 
* 


HS    SAY 


******************************************* 


BEBCFTf     ALLOCATION     A  KI    VAHIAELE       ASSIGNMENTS 


BEAL 
HEAL 
BEAL 
BEAL 
REAL 


BQE1J  (30C) 
-F{3CC] 

1CC,14) 


JITETF 
A 

„    PHMC1.2C1) 
REAL    TISC  (300) 
BEAL    LOFF.BU 
BEAL    IRRK1F  (30C) 
REAL    LI, 12 
BEAL    N ASCET,N0,ISCANG,ISOAN0, KA 

INTEGER    TITLE1  (60)  .TITLE2  (80)  ,TIIL£3  (80) 
IN1EGEB    EESCfi  (4C,E0) 

CCBBON    /AIBO/     EAIBO 

CCBMON    /EC/    TIFLK1,  VARBQ.  fcAVEEC,SCALEQ 
"i£3FLE/     EIA,EIAOES,  EEAtSZ,U0 


COBMON    /AE: 


CONS3AN1S     J.SSIGNMEN1S    ========================== =======C 


ai=20i 
jj=o 

TSTOLD=0. 
F4=1 

PI=2.*ABSIN  {1 . ) 
HLB=HL/2-H 
HAIHO=3. E4 
TAU=0. 

xns=i. E+e 

XCfl=1.  E+2 
IKS=1. E-3 
XK2=1.  E-7 
BID=90./AESIN  (1.) 
REARTH=6.4E6 
BU=3.986E14 

THETSF=2.*PI/ (2  4.*3  6  00. ) 
VSUBF=BEABTH*'IEE1£F 

CALL   ERESET     (208,256,-1,1,1) 

:=========    DATA    INEOT    SECTOB 


BEAD  (5,530) 
BEAD  (5,530) 
BEAD  (5,530) 


READ 
BEAD 

READ 
BEAD 
READ 
READ 
READ 


5,510) 
'5, 510 

5,510 
'5,510 

5,510 
|5,E10) 

5,510 


7I11E1  (I)  ,1  =  1  ,8C 
III1E2  (I)  ,I  =  1,8C 
TIIIE3  (I)  ,1=1,80 

(DESCR  (1,11)  ,11=  1,50) ,  CIA 

EESCR  (2,11)  ,11=1,50)  ,  DIACBS 

CESCE  3,11)  ,11=  1,50)  ,  EEAESZ 

EESCH(4, II  ,11=1,50),  V.AVE 

EESCR  5,11)  ,11=  1,50)  ,  ETCTAL 

i  EESCR  (6,11)  ,11=  1,50  ,  IHSEE 

(EESCR  (7,11)  ,11=  1,50) ,  RGRND 


(EITHER  TEFLMT  CE  HAVEBQ  BUST  EE  DEFINEE,  BUT  NOI  BOTH) 
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READ 
READ 
BEAD 
READ 
BEAD 
BEAD 
BEAD 
READ 
BEAD 
READ 
READ 
READ 
READ 
READ 
BEAD 


510) 

5  10) 

510 

510 

510 

5  10) 

510) 

510) 

510) 

510) 

510) 

510) 

5  10) 

520 

520 


[EESCR 
EESCB 
EESCR 
EESCR 
EESCR 
EESCE 
EESCR 
IEESCE 
[EESCB 
EESCR 
EESCR 
EESCR 
EESCR 
EESCR 
EESCR 


6,11)  ,11=  1,50)  , 
9.Il[,II=1<50[. 
10, II), 11=1, 50) 


11,11 
12,  II 
13,11 
14.11 
15,11 
16,11 
17,11 
18,11 
19,11 
20,11 
21,  II 
22,11 


,n 

,n 
,n 
,n 
,n 
,n 
,n 
,n 
,n 
,n 
,n 
,n 


1,50 

=1,50) 

=  1,50) 

=  1,50) 
=  1,50 
=  1,50) 
=  1,50) 
=  1,50 
=  1,50 
=  1,50 
=  1,50) 
=  1,50 


TDFIMT 
fiAVEBQ 

5CAL3Q 

HTRANS 

HSAT 

IHETHX 

LCFF 

RHCO 

VO 

SIGJIT 

AEAF 

AOELOtt 

AVGSPT 

NFIAGA 

NCISO 


II    FULL    AC    CCEEECTlCh   IS    HOT    ISED,    THE    following 
PABAMETEES    ARE     NCI     REALLY    USEE.    CNLY     BWIDTK 
HAS    AN     EIEECT.     it    CACSEs     the     ICCK    AHEAD    ANGLE    TO 
BE    LARGER    BY    TAU*EEIEOT. 


READ  (5,510)  (EESCR  (23,  II),  11=1,  50)  ,  ABSLOZ 

BEAD  (5,510)  (IESCB  (24,  II)  ,11  =  1,50)  ,  XJT 

READ  (5,510)  (EESCB  (25, II), 11=1,  50)  ,  B»IDTH 

READ  (5,510)  EESCR    26,11)  ,11  =  1,50)  ,  NA 


I3EBATICE    LOCi    IlfllS: 

N1    NUMBEE    OF    S0ETEE1A    INTERVAIS    FROM    0  IC    THETMX 

N2     NOHBEE    OF    IhTEGRATION     INTEEVALS    FOR  ABSORPTION 

N3    NOHBEE    OF     INTEGRATION     INTERNALS    FOR  TOREOLENCE 

N4    NOMEEE    OF     IMEGRATION     INTEEVALS    FCE  PTF 

N5    NOHBEE    OF     SUEINTEEVALS    OSEE    FOE    SLANT     PATH     OPDATE    FOR 

THERMAL    EIOOK1KG 


BEAD 
BEAD 
BEAD 
BEAD 
BEAD 


5, 520) 
5,5  20 
5,520 
5,  520 
":20 


l§: 


EESCR  (27, II) ,11=1,50)  ,  N1 

EESCR  (28,11) ,11  =  1,50    ,  N2 

EESCR  (29, II), 11=1, 50)  ,  N3 

EESCR    30, II), 11=1, 50    ,  N<4 

DESCE    31,11) ,11=1,50    ,  N5 


ECHC    CHECK    OUT  POT    SECTOR    ==  =  =  =  =  =  =  ===  =  ===========  =  ======c 


HBITE 
BBITE 
HBITE 
HBITE 
HRITE 
HRITE 
HBITE 
HBITE 
HBITE 
BBITE 
SBITE 
BBITE 
HBITE 
HBITE 
HBITE 
HBITE 
HRITE 
BBITE 
HRITE 
HRITE 
HRITE 
HRITE 
HBITE 
HRITE 
BBITE 
HBITE 
HRITE 
HRITE 


,190 
,  195 
,195 
,195 
,180 
,  180 
,180 
,  180 
,  180 
,  180 
,180 
,  180 
,180 
,  180 
,180 
,  180 
,  180 
,  180 
,180 
,  180 
,180 
,  180 
,180 
,130 
,  185, 
,185) 
,18  0', 
,180) 


(Till 

TI1L 

TITL 

EESC 

lEESC 

EESC 

EESC 

(EESC 

lEESC 

,DESC 

(EESC 

<EESC 

IDESC 

(DESC 

DESC 

DESC 

DESC 

(DESC 

(DESC 

DESC 

DESC 

(EESC 

(DESC 

DESC 

DESC 

DESC 

DESC 


"ill?' 

2,11 

3,11 
4,11 

5,11 
6,  II 

I'11 
8,  II 

9AJI 

10,1 

E  (11  ,11 
B  (12,11 
13,11 
14,11 
15,11 
16  ,11 
17,11 
18,11 
19,11 
20,11 

21  ,11 

22  ,11 

23  ,11 
24,11 


1,  EC 

I  ,80 
1,80 

II  =  1 
11  =  1 
11=1 
11=1 
11  =  1 
11=1 
11=1 
11=1 
11=1 
,11  = 
,11  = 
,11  = 
,11  = 
,11  = 
,11= 
,11  = 
,11  = 
,11  = 
,11  = 
,11  = 
,11  = 
,11  = 
,11  = 
,11  = 


»I2 

,50 
,50 
,50 
,50 
,50 
,50 
,50 
,50) 

1,50 
1,50 
1,50 
1,50 
1,50 
1,50 
1,50 
1,50 
1,50 
1,50 
1,50 
1  ,50 
1,50 
1,50 


EIA 

DIAOBS 
EEAHSZ 
HAVE 
PTOTAL 
1HSEE 
HGBND 
TDFLMT 
HAVEBQ 
SCALBQ 
HTRANS 
HSAT 
TEETHX 
LCFF 
RHCO 
VC 

SIGJIT 
ADAP 
ACBLOM 
AVGSPT 
NFLAGA 
NCI50 
ABSLOZ 
XJT 
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SEIT2  (6,180)  IDESCB  (25,11)  ,11  =  1,50)  ,  3SIDTH 

HEITE  6,130)  IDESCE    26,  II)  ,11=1,50)   ,  NA 

KRITE  6,185)  DESCH   27,11)  ,11=1,50)  ,  N1 

SRITE  (6,185)  DESCB  (28  ,  1 1)  ,  1 1  =  1 ,  50     ,  N2 

HEITE  (6,185)  (DESCH  (29  ,  II)  ,11=1 ,50)  ,  N3 

SRITE  6,185  (D E £C 3    30 , 1 1    , 1 1  =  1  , 50     ,  N4 

WBITE  (6,185)  DESCS    31  ,11)  ,11=1,50)   ,  N5 


C 

C==  =  ==========    DATA    REALIGNMENT    SEC1CB    ====  =  ===========  =  =====  =  =  =  ===  =  =  =C 

C 

If     (SCA1EC.2Q. C.)     SC AL3Q=CI A/5. 

II    JWAVEEQ.NE. C.)     V  A  ESi  =  (  2.  *  P  1*  Ji  A  VEB  £)  *  *2 

SGJIT0=SIGJI1 

If     (NFLAGA.  EC-  1)      AEAP  =  0- 

If     (NOISC.EQ-1)     I£CAK0=1.E10 

H1SQBT=SCHTTNA) 

L0GA3S=-AL0G (AESLCZ) 

TAU=1./EJiIDTH/EI 
C 
C=  =  =  =  =  =====  =  =  =    INITIAL    CALCULATIONS    SECTOB    =  =  =  =  =  ===  =  =  =  =  ========  =  =  =  =  =  =  =  C 

C 

If     (RHO0.EQ.  0.  ANC.THSEE.NE.  0.  )    BHO0=  ( 5iA  VE/ .  55Z -6)  **1.2*(.054/THSEE 

IHETHX=TEI1MX/ETD 

HSA1D=HSA1-HTEA*S 

HAEVGD=HTSANS-fcGEND 

IF     (HABVGD.LI.  C- )    HAEVGD=0- 
C 

R=REARTH+HTRANS 

RS=REAHTE+HSAT 

NFLGSN=0 
C 

C 

C  COMPUTE    EABTH    CENTER    ANGLE    OFI-SET 

C 

ANGOFF=LCFF/R 
A2=ANGOFf/2. 
C 

C  COMPUTE    COORDINATE     TRANSLATIONS    DUE    TO    CFF-SET 

C 

C 

Li  =  +  2.*R*£IN(A2)*COS(A2) 

1Z=2.*B*SIN (A2) **2 
C 

C 

C      COHPOTE  EARTH  CENTER  ANGULAR  RATE 

c 

C 

IBEIDT=SCBT(MD/B5**3) 
C 

c 

C      COMPUTE  CEEITAI  SFEEE 

C 

C 

V£AT=THETDT*RS 
C 

C  ADJUST    INTENTANECOS    SLEHRATE     fCB    SOTATICN    OF    EARTH, 

C  assuming    toth    are    cclinear. 

c 

C 

THETDT=AES (T HETIT-TEETS F) 
C 

c 

C  COMPUTE    VELOCITY    CF    SATELLITE    EELATIVE    TO    THE    TRANSMITTER    SITE 

c 

VS=THETD1*RS 
C 

C  COMPUTE     1NITIA1     LIMITS    ON    THE    TIME    AND     ANGLE    EARTH    CENTER 

C  .for    an    or   line    flight    path 

c 
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ECANG=ARCCS (R  *  SI  N  (TH  ETM  X  )  /R  S)  +1HETMX-PI/2. 
1IHIX2=ECANG/TE2TE1 
TIHTOI=2.*TIMEX2 
ETIHE=TI£EX2/N  1 
C 

C  If    SATELITE    ALITITUEE    IS     AS     TEE    OEDEE    Of    THE    RADIUS    OF 

C  THE   earth,    gutsff     (footprint)       WORKS    BETTER    BY    SETTING 

C  the    time    step    instead   ox    the    aigle   step 

c 

c 

IF     (HSAT.GT.  REABTE/2.)     CTIHE=10. 
C 

C 

C  CCHPDTE    VACUUM     IBEAEIANCE    AND     APERTUEE    MTF.        aCTE    THAT    THE 

C  PCINT    SPACING     EEEE    MCST     BE    THE    SAME     AS    USED    IN    THE    MTF 

C  INTEGRATION    3ELCW.  HENCE    WE    EEFINE     DX    NOW. 

C 

EX=DIA/N4 
C 

C 

C  THE    CALL    TO    D0CS1    DETERMINES     TEE    CONSTANT    THAT    MAKES    THE    EXIT 

C  AFEBTORE    EOHEB     ETCIAL. 

c 1 _ 

c 

CA1L  D0CS1  (PTCTAI,flJ.) 

CALL  FARfLD  (E X  ,  1ISC  ,IR R MTF , H I , NU ,HS ATD , W AVE , PTOT AL , PMAXO) 
C 

c 

C      IF  THE  1AVE3Q  IS  NOT  SPECIFIEE  AS  INPUT,  COMPOTE  ON  BASIS  OF 

C      TDfLMT. 

c 

IF     (WAVEEC..EQ.  C-)     CALL    DETWAV     (IERMT  F  ,  DX  ,  N4) 
C 

c 

C  IF    TDFLMT    HAS     fcCT    SPECIFIED,     COMPUTE     IT    . 

C 

c 


IF     (TDFLMT. NE.  C.)     GO    TO     10 

:t    BQIEEL     (IFEM' 
TEFLHT=SQET(  1./TECJ 


CALL    BQIEEL     ( I I i MT F , ^U, DX  ,TBQ ) 


10  CONTINUE 

C 

C  DEFINE    BEAM    QUALITY    MTF    ARRAY. 

c 

CALL    BM     (EQMTF  ,DX,N4) 
C 

C  DEFINE    JITTER     MTF    ARRAY. 

C 

c 

CALL    JIT     (JITMTf  ,  EX  ,  N4,  SG  JIT  C  ,  J.A  VE) 
C 

C 

C  COMPUTE    ZENITE    ABSOBETION,    SCATTERING,    AND    RHOO    VALUES. 

C 

C 

CALL    ABSCEB     ( K 2  ,HTB A NS. T AESO ) 

CALL    SCAT     (N2,HIEANS,TSCAT0) 

IF     (HHOO.EQ-0.)     CA1L    SHOTRB     (  N2  ,  HTEA  NS  ,  E  AE  VG  D,  H  AV  E,  RHOO) 
C 

c 

C      DEFINE  SEEING  AT  E50CA 

C 

c 

THSEE=  (iJ¥E/.£5E-6)**1.2*{.  05  4/BHOO) 
C 

IF     (NFLAGA.GE.  1.  AKD. NOISO.EQ-  C)     CALL     ISCTRB     (N 3, HTR ANS , HAB VGD , fc A VE 
1,ISdAN0) 
C 

C 

C      CALCULATE  LOG  AMPIITUDE  SCI  NT  1 1 LATIO K  FCR  ZENITH  ANGLE. 
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c 

c 

CALL    SCIM     (SIGXZ,N3,HTR  A  N3  ,  h"  GF  ND  ,  KA  V  E) 
C 

C — ■ 

C  CCBPOTE    TEE    TEIEB1I    ELOOMING    FfcASE    DISTCETION    THAT    does    not    depend 

C  SLANT    PATE    CH A E ACI EB 1ST ICS 

c 

C 

CA1L    PFTFES     (MI  ,M I B , 5 AVE , FH) 

CA11    HESEES       MI.KLH.FH) 

CALL    PHVAB     (MI,fiLH, F TOTAL ,P H  ,  SI GSQO) 
C 

c 

C  INITIALIZE    ACCUMULATED    FLUENCE    TC    ZEEC 

c 

C 

FLDX=0. 
C 

C 

C      SET  SNB  FLAG  TC  CFF 

c— 

C 


NFLGSN=0 

C 

c 

r ***»*»»**** ** *     MAIM  PHCGFAM  LCOP     ***************** Q 

C 

c 

C  LCOP    ON     IEEADIANCES     AND     ACCUMULATED     FLUENCES    FOE    A    FIXED    LOFF 

DO    120    1=1,  N1 
C 

c 

C  COMPUTE    TCTAL    LINGTE    OF    ILLUMINATION    TIKE    F30M     (+,-) 

C  OVEB    PASS    ANGLE 

C  THIS    IS     THE    TCTAL    TIME.        SE     USE    HALF     OF    THIS    TO    DETERMINE    THE 

C  FUNCTION     EVALDATICK    ANGLE. 

C 

c 

TIME=2.*  (  (1-1)  *ETIMI-»DTIME/2.  ) 
C 

C 

C  COMPUTE    EARTH     CENTEE    ANGLE    AT    CN    LINE    CCOBDINATES     AND    TIME/2. 

C  THIS    ANGLE    IS     TC    TEE    MID-POINT    OF    THE    INTEGRATION    INTERVAL. 

C  THE    ANGLE    £C A N GF, CCM FUTED    LATEE,    IS     TO    THE    UPPER    LIMIT    OF    THE 

C  integration    liHit. 

C 

C******    1HE    COMMENTED    ECANG    SHOULD     EI    USED     WREN    IT    IS    DESIRED    THAT 

C  THE    INTEGEATICN    SHCULD     STAET    AT    THE    LCS     POINT    AND    INCREASE 

C  TOWARDS    THE     ZENITH.       THE    CUEFENT    CALCULATION    STARTS     AT 

C  THE    ZENITH     A  K  E    CCES    DOWN. 

C 

C***   ECANG=  (T1BT0T-TIME)/2.*THETDT 

C 

ECANG=IIFE/2. *TEETDT 
C 

C 

C  COMPUTE    TARGET     COCEDINATES 

C — • 

C 

X0=BS*SIN  (ECANG) 

Z0  =  RS*COS    ECANG) 

X  =  X0 

Y  =  IY*COS  (ANGCFFJ  +  (ZO-B+LZ)  *SIN|ANGOFF) 

Z=-LY*SIk  (ANGCFF)  +  (ZO-E  +  LZ)  *CCS  (ANGOFF) 
C 

C 

C  COMPUTE    CN    LINi    OF    SIGHT    ANGLE    OF    SAT. 

C  THIS    ANGLE    IS     NCT    USED    AT    PSESENT    .  IT    IS    THE    ANGLE    CF    THE 

C  SAT    AS     MEASURED    FECM    A    POINT     ONDEE    THE    GROUND    TRACK     . 

C 

c 

IHETA=ATAH  (RS*SIN(ECANG)/(RS*CCS(ECANG)-R)) 
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c 

c 

C      COMPUTE  RANGE  1C  1AEGET 

c 

C 

BANGE=SQ£T (X**2+Y**2+Z**2) 
C 

C  INSTANTANEOUS     £IIS    BATE 

c 

C 

VX=VS*COS  (ECANG) 

VY=-VS*SIN  (EC  A  KG)  *SIN  (ANGCFF) 

VZ=-VS*SIN (ECANG)  *CCS(AHGCFFJ 

IX=Y*VZ-2*VY 

HY=Z*VX-X*VZ 

WZ=X*VY- Y*VX 

PHIDOT  =  SC.BT (KX**2+WY**2+MZ**2)/BANGE**2 


C 


C  CCHPOTE    ANGLE     IKE    ZENITH 

C 

C  OMEGA     EQUALS     ANCIE    DCWN     FROM     ZENITH 

C  CMfcIND     EC.0ALS     AKGIE    CF    WIND     ATTACK    TC    LCS    I?    TARGET 

C  MOTION     ANE    WINE     ABE     COPLANAB. 


OMEGA= AB CCS (2/ BANG E) 
CHHIND  =  ABSIN  (X/RANGE) 


COSOMG=CCS(OMEGi) 
SECOMG=1 ./COSCPG 
CCSBND=CCS(OMSIND) 
C 

C  ADJUST     FCE    SLANT     IATH,     BLOOMING    LOSS,     AESCBPTION     LOSS,     SCATTERING 

C  LOSS,     ANE    TDBEDIENCE    RHO 

C 

c 

C 

CALL    AV      |H5#VC,EHICCT,HTBANS,CCSHND,COSOBq,E) 
SIGSQ=SIGSQO*E 

c 

C      APPLY  THEBMAI  ELCCK1NG 

C       ADAPTIVE  OPTICS  DEGREE  OF  COMEENSATION 

C ■  — 

C 

SIGSQ=SIGSQ*ACE1CE 
C 

SIG=SQBT  (SIGSC) 

TELOOM=EXEi-SIGSQ) 

IF     (SIGSC..C-T.  1.2)     CAIL    3L0OM      (EI  A,  BE  AMS  Z  ,  SIGSQ  ,T  BLOOM) 

TAES=TABSO**SECCMG 

TSCAT=TSCATO**£ECCHG 

BHC=BHO0*COSCMG**  (.6) 
C 

C 

C  COMPUTE  LOOK  AHEAD  ANGLE  ASSUMING  EAETH  BOTATION  AND  SATELLITE  MOTION 
C  ABE  IN  THE  SAME  DIRECTION. 

c 

c 


TBCT=RANGE*2./3.EE 

DBI=THCT*  (ABS  (^SAl-VSURF)  ) 

PBCJAN=AECOS  (X/EAKGE) 

PBCJ=DBT*SIN    EECJAN) 

ND=PROJ/EANGE+TAU*PHIDOT 

IF     (I.EQ.1)     ZENNU=NU 
C 
c 

C  IF   executed,     the    following    is    a    full    ao    simulation 

c 

C 


IF     (NFLAGA.EQ. C)     GC    TC    40 
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c 

C      HCDIFY  EEC  FCE  EIIECT  OF  ADAPTIVE  OPTICS 


C 


C 


BHCLD=BHC 
RBOU=BHO 
E0=2.  1  *BHC 


C  DE1EBMINE    APFECXIEA1E    I-EEL    TC    SEE    If     AC    SHOULD    BE    USED, 

C  given    PEESENt    NOISE. 


DCE=DIA/EHO 

CALL    TBBEEL     (I E E Mil , EIA , N H, D X  ,  EXDTST ,SHC) 


C  COMPUTE     BESIDUAI     VARIANCE    DOE    TC    PEEFECT     ADAPTIVE    OPTICS, 

C  INFINITE     EANDVilETfc,     AND     FINITE    NUMEEE    OF    ACTUATORS. 

c 

E1  =  -32  0*  JEIA/EO/NASQET)  **  (1.6  6  67) 

TES=LOGAES*SECCHG 
TSENSR=E1F  (-IE£) 

c 

C  DETERMINE    PHASE    VARIANCE     ASSOCIATED     SITE    SENSOB. 

PHSERB=8.5E-6*EKIETH*NA**2* (RANGE/5.  E5*2. 5/DIA) **4*  (  80 ./I  J  T* . 5/1 SE 
1NSE)  **2 
c 
c 

C  ISCBEASE    EESIEDAI    PEASE    VABIAbCE     AND     CCPPUTE    THE    I-EEL. 

D1=D1+PHSERB 
C 

C  COHPUTE    I-EEL    CI    AC    TUBBOLENCI    CCREECTED    BEAM. 

C 


C 


EXD=EXP  (-D1) 


C  IF    AO    COEBECTEE    I-EE1    IS    LESS    THAN    NC    AC    CORBECTION,    ASSUME 

C  ONLY    TILT    COHEFC1ICN    IS     USED. 


If     (EXD-I1.EXD1S1-OB.NFLGSN.EC-1)     GO    TO    20 
CALL    FNDEBO     {  E XE  ,  EHC , DI A , EHOO  ,£X , N4 , IRE HTF) 


GC    TO    30 
20  CONTINUE 

C 

c 

C  THIS    SETS    A    F1AG    TC    FEINT    THAT    AO    SYSTEM    IS    ONLY    TILT 

C 

RBC=RHOLE 

HI IGSN=1 
C 

30  CONTINUE 

C 

C ■ 

C  ESTABLISH    TEST    INCBESENTS    TO    SEE    IF    FULL    ISOPLANATIC    CALCULATION 

C  NEEDS    TO    EE    DCNE. 

C — 

C      COMPUTE  TEE  ISCELANATIC  ANGLE. 

c 

C 

ISCANG=CCSOMG**  (2.6667)  *I£OAN0 

TSINES=NC/ISCANd 

TSTDIF=AEE  (TS1NEV(**1  .66  66  67 -T  ST  OLD) 
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ETSDIP=EXF (TSTEIF) 

IF     {NISO.EQ. 1)     ISIOLE=TSTNEW**1. 666667 
c 
c 

C  EETEBMINE    IF    ISCELAKAHC    MTF    SBCULD    EE    BE-COMPUTED 

c 

c 


NISO=0 

iISO=1 


IF     (ETSDIF.GT.  1.  1)     N: 

II     (I.EQ.l)     NISO=1 

if    (noisc-eq.1)    niso=0 
40  CONTINUE 

C 

C 

C  SEDUCE    PEAK    IEB1D3AKCI    BY    SLAM    DISTANCE 

c 

C 

IMAXSC= (ESATE/BANGE) **2 
C 

C    IBIS  ALLOWS  ONE  TC  INCLUDE  TILT  ISOPLA N AT  IS B.   NOTE  HOWEVER,  THAT. 
C    IF  THE  FULL  ADAFT1VI  OPTICS  ALGORITHM  IS  CSED  WITH  THE  ISO? L A N ATIS K 
C    CCEEECTION,  THEN  THIS  CALL  SHCUIE  NOT  EE  US  ED- REPEAT- SHOULD  NOT  BE 
C    USED.    ALSC  NOTE  TEAT  A  VEETICAI  1URB0LENCE  PROFILE  IS  NEEDED 
C    SiHEN  EVEB  THIS  FCOTINE  IS  USED. 

c ' 

c 

IF     (NFLGSN.EC-  1. AKD.NOISO.EQ.O)    GO    TO    5C 

IF     (NOISC.EQ.  1.CB.NFLAGA.  EQ.  1)     GO    TO     60 
50  CONTINUE 

OMX=NU 

CMY=0. 

CALL    TLTISO     { E I  A , CMX ,03 T , BESTIT , SECO HG, HTE A3S , HAB VGD) 
C 

SIGJIT=SCBT (SGJITC**2+aESTXT) 
60  CONTINUE 

C 

C  COMPUTE     IEBADIANCJ     EECEEASE     DCE    TO    AMPLITUDE    SCINTILLATION     FOE    THE 

C  OFF    ZENITH    PASS-       NOTE   THAT    TE1S    AS    LOSS    HILL    ONLY     EE    APPLIED    IF 

C  FULL    ADAPTIVE    OPTICS    COMPENSATION    IS    ASSUMED. 

c 

CALL    SINT1S     (TAMP,SIGXZ,SECOHG) 
C 

C  CALCULATE    EFFECTS    CF    JITTEB     AN!    TUEBULENCE 

C 

S0fl=0. 
DX2=DX/2. 
C 

C  IT    IS    IMECKTANKT     TO     NOTE     THAT     ALL    THE    ASSAYS     ARE    DEFINED 

C  TO    BEGIN    AT    DI/2. 

c 

X=DX2 

DO    110    J=1,N4 

IF     (X.GT-EIA)     GO    TC    100 

F1=IBBMTI  (J) 

CALL    MTFATM     (X  ,EI  A  ,  BHO,  AD  AP  ,  F  2) 

F3=JITMTF  (J) 

IF     (NFLAGA.E^.  LAND.  NFLGSN.  EQ.O)     GO    TO    70 

IF    JNOISC-EQ.O)     CALL    MTFJIT     (  X  ,SIG  JI I  ,  S  A  VE  ,  F3) 
70  CONTINUE 

FU=1- 

IF    JNFLAGA.EQ.G.CE.  tiCISO.  EQ.  1  .CE- NFLGSN.EQ.  1)     GO    TO     90 

F4=TISG (J) 

IF     (NISO.EQ.l)     CALL    ISOPLA     ( J J  ,  NU , X, SECCMG , R AV E, FU , HTB ANS , H AB VGD  ,T 
1ISC.NU , J) 
90  CONTINUE 

F5=BQMTFJJ) 

SUM=SUB+F1*F2»F3*FU*E5*X 
100  CONTINUE 

X=X+DX 
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110  CONTINUE 

TTUBB=S0B*2. *FI*DX 
C 

C  CALCULATE    THE    EIS011IHG     I NSTS NT AN EOU £    INTENSITY 

C  CCNVEBT    TC    AN     APFECCE    WHICH     AVERAGES    THE    RESULT    OF    AN    ESS 

C  TBEATMENT    OF    THEBCAI    BLOOMING    WITH    A     MULTIPLICATIVE    TREATMENT 

C 

c 

SSC.T=1./TTURE-  1. 

S£QB=1 ./IELOCK-1. 

SSC=SSQT+SSQE 

TKULT=TTUE3*TE1CC2 

TBSS=1  ./ (1.+SSC) 

CALL    BELIC1     ( 1 1CT AI , TESS , 1MU LT ) 

PMAX=PMAXO*TAES*l£CAT*TMAXSC*'nCTAL 
C 

C  APPLY     SVEEAGING     EACTOE 

C  APFEOPEIATE    TC    SPC1     £I2E     CESIEED 

C 

c 

PSAX=PMAX*AVGSET 
C 
C 

C 

C  IF    FOLL    AC    IS     USED,     EECEEASE    IEEADIA NCE    EOS    AMPLITUDE 

C  SCINTILLATION     EFFECTS 

C 

c 

TAHPL=  1. 0 

IF     (NFLAGA. EQ.  1.AKD.  NFLGSN.  EQ.O)     TAMPL=TAMP 

IF     (TAMPL.LT.  0-5)     TAEPL  =  0.5 

FHAX=PMAX*TAMEL 
C 

C  ACCUMULATED    FDII 

C 

c 


FI0X=FLU2*PMAX*CTIflE 
C 

c 

C  COMPUTE    IHTSGEAIED    IIM2S     AND     ACCUMULATED     ANGLES.     THE    IREADIANCE 

C  FUNCTION    HAS    EEEN    EVALUATED    AT    INTERVAL    MID-POINTS. 

C 

c 


ECANGF= (1IME+E1IME)  /2.*THETDT 

XOr=BS*SIN(ECAKGF) 

Z0F=BS*CCS (ECANGF) 

XJ=X0F 

YF=LY*COS 1ANGCFF 

ZF  =  -LY*SIF  (ANGC~ 

BANGEF=SCET  (Xf  " 


FF)  +{Z0F-R  +  LZ}* SIN (ANGCFF) 
CFF) + (20F-B+LZ) *CCS (ANGOEF) 
**2+YF**2+ZF**2) 


OHEGAF  =  AECOS  (ZI/BANGEF) 
C_ 

C  TIME    OF    IILDMIliATICK    FBOB    ZENI1H    TO    TEETA.     THIS    IS    ALSO 

C  EVALUATED    AT    TEE    IOHEB    LIMIT     JND    NOT    THE    MID-POINT. 

C  IF    THE    COMPUTATIONS    ARE    DONE    EEOM    THE    THETMX    TO    ZENITH    INSTEAD 

C  OF    THE    WAY    THEY    ABE,    THE    TIME    DEFINITION    NEEDS    TO    DE    CHANGED. 

c 

IILLUM=TIME/2. +ET1ME/2. 
C 
C============  =    LOOP    BESUL1S    STCBAGI    SECTOE    ================ ===========c 

C 

A  (I,1)=I 

A  (1,2)  =BABGE*XKB 

A    1,3    =PEIDOI/ XKW 

A    1,4) =OKEGAr ♦BID 

A    1,5) =TIILUM 

A     1,6    =TECAT 

A  jl,7)  =TAES 

A  (1,8)  =EEC*XCR 
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A  (1,9)  =TT0R3 

A  (I, 10) =TAM?L 

i  1,11  =IELOCB 

A (1,12) =lflAXSC. 

A  I, 13)=EKAX*XK2 

A (1,14) =FLUX*XK2 
120    CONTINUE 
C 
0==============  FINAL  EESUITS  OUTP01  SECTOR  ===================== ======C 

C 

E  (1)=DIA*XCM 

E  (2)=WAVI*XMW 

E  (3  =PIC1AL/XBS 

B  (4)=TDFIBT 

E  (5)=HTRANS 

E  (6)  =  HSA'1*XKW 

E  7  =THETBX*BTE 

B  8)=LOFF*XKW 

E  9)  =RHO0*XCM 

B  1U)  =V0 

E  11)=BGEND 

B  12)  =  IHSEE 


HEITE  (6,220) 

KBITS  6,240)  <  (A  |I,J)  ,  J=1,8)  ,I=1,N1  ) 

WRITE  6,230) 

WBITE  (6,250)  ( A  (I,  1  ) .  I A  (I . J)  ,  0  =  9  , 1 4 )  , 1=  1 , N 1 ) 

HBITE  6^330)  (Ejll)  ,U.=  1 ,  12) 

BEITE  (6,130)  ACE1CB,AVGSPT 

EIA=DIA*XCM 

DIA03S=DIAOBS*XCH 

EEAMSZ=EEAMSZ*XCM 

WRITE     (6,140)      EIA,EIA03S, EEABS2 

SCALBQ=SCA1BQ*XCE 

HRITE     (6,150)     SAVIEQ,SCALEQ 

B  (1)=ABSLCZ 

B  (2)=XJT 

B  (3)  =BHIETH 

B  (4)=NA 

WRITE     (6,380)      (E  (LI)  ,LL=  1  ,4) 

WRITE       6'260l 

II     (NFLAGA.EQ.  CAKE.  NOISO.EQ.  1)     WRITE     (6,310) 

IF     (NFLAGA.EQ.  LAKE.  NOISO.EQ.  0)     WRITE     (6,280 

IF     (NFLAGA.EQ.  LANE. NOISO.EQ.  1)     WRITE       6,290 

II     (NFLAGA.EQ.  0. AKE.  NOISO.EQ.  0)     WRITE       6,300 

WHITE     (6,340)      S  1  ,  }>  2  ,  N3,  N4  ,N5 

ESIGXZ=EXE f-SIGXZ) 

WRITE     (6,160)      SiGXZ,ESIGXZ 

ZENNU=Z£NNO*XMS 
WBITE     (6,320)      ZENK0 

RESTLT=SCET(RI£ILT) *  XMW 
WHITE     (6,360)     EESTL1 

SGJITO=SGJITO*XMS 

WBITE     (6,350)      SGJITC,ADAP 


FL0X=FLUX*XK2*2. 

TILLUM=TI1L0M*2. 

WHITE  (6,370)  FLU2.IILLCM 

II  (NFLGSN.EQ.  1)  WRITE  (6,170) 

S20P 
C 
C====  =========  I/C  FORMAT  STATEMENTS  ====  =========  =  =  =========  =  =  =  =====C 

C 

500    FCRMAT  (50X.F20.0) 

510    FORMAT   50A1,E20.0) 

520    FORMAT  (50ALI20) 

530    FORMAT  (6CA1) 

C 

130    FORMAT (1 CX,' BLCCEING  ADAPTIVE  CPTICS  FACTOR   =  ',£10.4,/, 


103 


140 

150 

160 

170 

180 
185 
190 
195 
220 


230 


240 

250 

260 
270 

280 

290 
300 
310 


320 

330 


340 


350 


360 

370 


380 


|  1 CX, ' IEBADIANCE    ABEA     AVIEAGING    FACTO* 

FCBMAT  (10X, 'TELESCOPE    DIM  EN  SIGNS:',/, 
J  1CX, 'CUTER    EI AMETEE 

1CX, 'INNER    EIAMETEE  = 

10X, ' GACSSIAN     WAIST     BIABETER    = 
FORMAT  (1 CX,  ' EMS     SAVES    DISTORTICN 
|  1CX, 'PHASE    CCREELTATION    LENGTH    -        ^io.h. 

FORMAT     (10X, 'ZENITH    LOG    AMPLITUDE    VARIANCE    =',212.4 
1  'TIVE     IRBADIANCE     5EDDCTICN    =  '     E12.4,/) 

FORMAT     (  1QX,  '  NCTE==  =  — —THE    AEAP7IVE    OPTICS    SY 
7IEE    DEGRADED    PERPCSMANCf    OVER    TEE 


S10.4,/) 


,F10.2,'     CM',/, 
,F10.2,  '    CM'  ,/, 
,F10    2,'     CM',W 
'  .E10. 4,/, 

E10.4,'    CM',/) 

/,10X, 


REIA  ' 


1QI,50i1, E20.O) 

Fill:' ,//) 
ANALYSIS     RESULTS:*,//, 


'ILT 


ITEM 
ONLY 


USED 
CASE' 


WILL 
,/) 


1CX,5CA1,I2Q) 
1E1,9X, 'INPUT    DATA 


in  i  ,  7  a  ,  ■ 

CX,80Al,/) 
1  H  1,91  -  'PATH 
6X,'STEE' , 10X 
5X,'ATMCS.     1 
7X,  'NO.  '  .  1  1X 


i    i  *  .       ia  .1 )    '      D 

6X,*  (SEC)  '  ,8X  ,'3CA7f  , 
FORMAT  (1E1,9X,  'THAN"" 


a     ANaLIilb      riiUUi:',//, 

OX    'RAN GE ', 9X, ' SL Ei*  ,6X, '0 MEGA ',7X,' TIME' 

TRANSMISSION'  ,  10X,  'RHC  ,/, 

X     •  (KM)   '     5X  ,  MMRAD/3EC)  '  ,4X,  '   (DSG)  '  , 


v  ,  '  (MBAD/3EC)  • 
SX^ABSOEP' ,9X 


(CM) 
-MISSION 'ANALYSIS    OUTPUT:  *  ,//k 
16X\  'TBANSBISSICN    CO  EFF ICI 1  NTS  '  ,  1  3X  ,  '  R  ANGZ  '  ,  SX,  '  M  A  X    IF.RAD' 
• FLUENCE' ,/, 

IT",  ft 


,    'II 
1   II. 

4X     »' 


DJ.'^'ti'.ii/i/'JE  TO',5X,'AMP  LCS5',4X, 
■  TEERMAI  '  ,4X  'SCALE'  ,5X, '  (KM/CM  2)  '  ,5X, 
•  (KJ/CM2)  ',/,  lOx,  '  SPREAEINd'  ,  6X.  'WITH  A 


•  (KJ/cn2)  ' ',/ .  Sox",  •  spr£aeincV'  ,  ex.  'with   ao*  . 5x,  'blooming'  ,//) 

OHM  AT  (7X,?3.C,5X,F1C.3.5X,F3.4,5X,?6.2,5X,?6.2,oX, 

F7.545X,F7.5.5X,F8.<I,/) 
OEM  AT  (7X,F3.C,5X,F10.3,5X,Flt.5,5X,F10.5,5X,F10.4,5X,E10.4, 

5X.E1G. 4./) 


FORMAT 
FORMAT 
FORMAT 

I 
FORMAT 

! 
FORMAT 

I 
FCBMAT 


i6x, ' TEIS  EON 

12X, ' INCLOEIS 

12X, 'INCLOEIS 

•ISOPLANATISfi 
(  121,  ' INCLODES 

'ISOPLANATISfi 
(  12X,  'INCLUDES  FULL 

•ISOPLANATISP") 
(121 ,'IS  TEE  EASIC  CODE  WITH  CNLY  SOME  M 
TILT  CCBEECTION' ,/,12X, 'AND  WITHOUT  AN 


THE  STATISTICAI  CICUL  MODEL') 
THE  FULI  2CNAL  AC  MODEL  WITH  AN 
MODEL1  ) 

THE  FULI  ZONAL  AC  MODEL  WITHOUT 
MODEL*  ) 

TII7  CORRECTION  WITH  TILT  ' 


AN 


T  (10 
,T  ( 1  E  1 


OX 


FOBMAT 

FCRMA_ 

1CX 
10X 
1  CX 
1CX 
10X 
1CX 
1CX 
1CX 
10X 
10X 
1CX 
1CX 
10X 

FOHMAT  (1 CX 
12X 
12X 
12X 
121 
121 

FCBHAT  (1CX 
10X 
10X 

FCBMAT  (1 CX 

FORMAT  /// 
1CX 
1  C  X 

FOBMAT  (10X 
12X 
12X 
12X 
12X 

ENE 


K0DIY«KA 


iASURE    OF     '  , 
ISOPLANATIC 


,'iinxn    LOOK-AHEAI    ANGLE 
9X,'CASE     LATA     AND    CALCULATED 
TBANSEITIER     DIAMETER 
CAVITY     SAVELENGTH 
APERTURE    TOTAL    POSES 
TIRES    DIFFRACTION    IIMIT       = 
(WIDE    ANGLE     SCATTERING)  •   ,/, 
7BAKSEIITER    ALTITUIE 
SATELLITE    ALTITUDE 
MAX     ZENITH    ANGLE 
FLIGHT    PATH    OFFSET  = 

E  H  C  0 

WIND    VELOCITY  =     ', 

HEIGHT    OF    GROUND  =     ', 

OPT    SEEING    AT    5500A  =     \ 

NUBEEE    OF    INTERVAI    STEPS     FOR: 
ANGULAR     INTERVAL 
AEECEET2CN    INTEGRATION 
RHC    CALCULATION 
MTF    CALCULATION 
THERMAL     ELOOHING 
2-S1G3A-E    BEAM    JITTER 
TURBULENCE    JITTER    REJECTION 
(BISIEUAL    =     INITAI    *     ADA  F)  '  ,/) 
ISOPLANATIC    JITTER     (2-SIGMA-P)     = 
OX. '=  =  "===  =     PROPAGATION     RESULTS 
INTEGRATED    FLUX    ON    TARGET 
TCTAL    ILLUMINATION    TIME 
FOB    TEE     FULL     ZONAL    A-0    MODEL: 
ZEKITE    TRANS.     AT     t-C    SENSOR 
TABGET    RADIANT    INTENSITY 
ADAPTIVE    OPTICS     BANDWIDTH 
NOfEEF    CF    ACTUATORS 


=     '  ,E10.  4, '     UBAD'  ,/) 
FACTORS: ' ,//, 
'  ,  F10.2,  '    CM' ,/, 
',F10.6,  '     UM  ',/, 
',F10.4,'     MW*,/, 
',F10.2,/, 

*,F10.2,'  M*    /, 

*,F10.2,«  KM*./, 

',F10.2,'  DEG',/, 

*,F10.2,'  KM',/, 

•,E10.4,«  CM't/. 

»    F10.2,'  M/SZC',/, 

',F10.2,'  M',/, 

E10.2, '  ARC-SEC ,/) 

''='•,110,/, 

=    ',110,/, 

=  ; ,iio,/, 

=    ',110,/j 

=    *,F10.2,'     URAD',/, 
" ,E10.4,/, 


,E10.4    /) 

^15.6,*    KJ/CH2',/, 
=    ',F15.2,'    SEC*) 

{,E10.4,/, 

',£10.4,'     W/STER' ,/, 
',E10.4,'     HERTZ',/, 
'  ,F10.0,/) 
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SUEROUTINE    ALFA     (1,1) 
C 

C CALCULATE    THE    1CTAL    BOLECOLAB    AESCEPIICN 

C 

DIBENSIOli    ALT  (20)  , ATA  (20) 
C 

£»**♦**  THIS    IATA     FCR    CO    5?  9     FHCH     AFKL    INFO,     REPORTEDLY     EASED    ON     AN 

C******  UPDATE    OF     TEE    SCCLATCHY     LINE    DATA.       THE    WAVELENGTH    FOE    THIS 

C******         TRANSITION    IS    4.99210    HICRCNS. 
C******  ATBOSPHERE     IS    HI E-L ATI! OE E    SUMMER,     CLEAR    DAY. 

DATA    ALT/0.0.  1.0,2.0.3. 0 .a. 0,5. 0,6. 0,7. 0.6. 0,9. 0,10. 0,12. 0,14.0, 
1  16- 0  -  18. C. 2 0.0 ,25. 0,30.0- 35.0.40.0/ 

DATA    ATA/C. 356,0.269  .0.  146, 7. 28E-02, 3.3  5E-02,1 . 5 1E- 02, 6 . 33E-0 3 , 

1  3.40E-C3,1.61E-03,  1.79E-03,3.46E-04,6.5OE-04,  1.78E-0  5, 

2  1.06E-05,6.21E-06,3.61E-0e,C.0,G.0,0.0,0.0/ 
NI=19 

IHDX=1 

A  =  C. 

IF     (  (Z*3.279)  .Gl.  1.E5)     EETURN 

H=Z*1. E-3 
C 

C SEARCH    FCR    THE    AIIITUDE    INDEX 

C 

DO     10    1=1, NL 

IF     (H. LT.ALT (I))     GO    TO    10 

INDX=I 
10  CONTINUE 

R=  (H-ALT  (INDX)  )/  (ALT  [INDX  +  1 )  -  AITJIMD  X)  ) 

A  =  R*  (ATA  (lNDX  +  1)  -ATA  (INDX))  +ATA  (INDX) 

A=A/T000. 

EETURN 

EBC 
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c 
c 

c 

c 
c*** 

c  *** 

c  *** 


SDEROUTINI    A1?S      (S,2) 
CALCULATE    THE    1CTAL    SCATTERING 


DIMENSION    ALT  (20)  , ATS  (20) 


**  *  TH 
***  OP 
***  TR 
***  AT 
CAIA  A 
1 

DATA    A 
1 
2 
3 
NI=19 
S  =  C. 

If    iiz 

E=Z*T. 
INDX=1 


IS    DATA 

DATE    OF 

AN  SIT  10 

MCSPHER 

LT/0.0  . 

16-  0 

TS/1.U  1 

2.60 

2.62 

3.28 


*3.279)  .CI.  1.E5)     RETURN 
£-3 


C 

C SEARCH    FCE    THE    ALTITUDE    INDEX 

C 

DC    10    1=1, NL 

IF     (H. IT. ALT  (I) )     GO    TO    10 

INDX=I 
10  CONTINOE 

R=  (H-ALT  (INDX)  ) /  (ALT 

S  =  c*  [ATS    INDX  +  1)  -ATS 

S=S/1000. 

RETURN 

EKE 


INDX  +  1)  -AIT  (INDX)  ) 
INDX))  +ATS  (INDX) 
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SUBROUTINE  A3SCE3  (N,HT,T) 
CALCULATES  ZENITH  THE  TOTAL  INTEGRATED  MOLECULAR  TRANSMISSION 


10 


COflflON    /kino/     EAISC 
HAT30-HT)/N 
=HT+DELE/2. 


DEIH= ( 


HEIGH 

ALPZ=0 

DC     10 

CALL    A 

ALPZ=A 

HEIGHT 

CONTIN 

A1EZ=A 

T  =  EXP  ( 

EEIURN 

END 


1=  1  ,N 

LEA     (ALF,fcElGET) 

LP2+ALP 

=HEJGHT+EELE 

UE 

LPZ*DELE 

-ALEZ) 
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SCEEO0TINE  SCAT  (K,KI,T) 
C 

C CALCULATES  ZENITH  THE  TOTAL  INTEGRATED  SCATTERED  TRANSMISSION 

C 

COMMON    /A1MO/     KATtC 

DILH= (HAIMO-HT)/N 

HEIGHT=H3+DZIB/2. 

AIPS=0. 

DO     10    1=1, N 

CALL    ALFS     (ALS,EEIGH1) 

AIES=ALPS+ALS 

EEIGHT=HEIGHT+£ELE 
10  CONTINUE 

ALPS=AL?S*DELE 

I=EXP(-ALES) 

RETURN 

END 
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c 

c— 

c 


SUEfiOUTINE    EHC1FE     (  K  ,HT , HG, W , EEC) 
COMPUTES     XUF.AS     ZEKI1E    ATMCSPKEEIC    COEEEENCE    DIAMETER 


10 


CCHMON 

TPI=4. 

CK=IPI 

CKSQ=C 

DEIH=j 

HEIGHT 

SDH=0. 

DC    10 

CAIL    C 

S0B=SU 

HEIGHT 

CONTIH 

CM  =  SU 

CEN=2. 

EHC=]6 

EHC-RH 

EETOBN 

ESE 


/A1SO/     fiAlCC 
♦AESIN  (  1.) 

K*CK 

HAlflO-KT) /FICA1  (N) 

=HG+DELE/2. 

1=1, N 

N2E    {HEIGHT, CH2) 

B+CM2 

=HEIGH1+EEIE 
OE 

B* 
91 


0/ 


EELH 

*CKSQ*CNL 
8/DEN)  **  (-6) 
2*  1 
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10 

c 


SUEROUTINE    FAEFLD     (C 12, F ,  G, M I , lU , HSA ID , «A V Z, PTOTAL, PHA XO) 

COMPUTES    THE    FAEFIEir    IBRADIANCE    0?    TEE    APERTURE    DISTRIBUTION 
ADD    THEN    COMPUTES    TEE    COF.  RESPC  NDI  NG     APERTURE    Htf 

DIMENSION    F  (MI)  ,G  (Mil 

CCMMON    /ARBFLD/    LI I , tlAC ES, BE  A ESZ ,U0 

EI2=ARSIN  (1.) 

Pl=2. *PI2 

TSCPI=2. *FI 

CK=TWOPI/SAVE 

CKB=CK/HSATD 

DIASQ=DIA**2 

DCESSQ=D110BS**2 

R FOCUS =1 . 2*HSATD*SAVE/3QET ( DI 2 Sy-DOB SSQ) 

FACT0B=4. 

FACTOR=2.5 

FACT0R=5. 

RIIH=FACTCR*RFCCCS 

DXO= (DIA-EIA0ES)/2./KL 

DX1=ELIH/BI 

E1=DX1/2. 

DC    20    1=1, ML 

SUM=0. 

B0=DIAOBS/2. +EX0/2. 

DC    10    J=1ZML 

Z=CKR*B0*E1 

CALL    JO     (2, A) 

CALL    FIELD     (BO. OF) 

SDB=Sna+A*UR*BC 

E 0=  R  0+  DX  C 

BY  DIVIDING  EY  ETCTAL,  THE  FTF  Will  EE  UNITY 
F  m=  (SUM*CKF*DXO)  **2/P TOTAL 
£1  =  B1+DX1 


AT  THE  ORIGIN 


20 
C 

C COMPUTE  EEAXO  CK  AXIS 

C 

SUM=0. 

B0=DIAOBS/2. +DX0/2. 

DC    30     1=1, ML 

CALL    FIELD     (BO, OH) 

SOM=SUM+GB*RO 

EO=RO+DXC 

PMAXO=  (SUR*CKE*DXO) **2 


30 

C 
C— 
C  — 
c 

C 
C  — 

c 


40 
50 


CNORM-     PIXES     TEE     ETF    SO     THAT     TEE    MTF     INTEGRATION     PRODUCES    AN 
IBEL    VALUE. 

CNOBM=  (VAVE*BSATD)**2*PMAX0/PTCTAL 

COMPUTE    TEE    KTE 

R2=DX2/2. 

DC    50    1=1, MU 

SUM=0. 

B1=DX1/2. 

DC    40    J= 1.ML 

Z=CKB*R1*E2 

CALL    JO      (Z.A) 

SUM=SUM+A*F  (J)  *E1 

E1=R1+DX  1 

G  fl)  =SUM*TWOPI*DX1/CNORB 

R2=R2+DX2 

RETURN 

END 
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SOEEOUTINE  JO   <X,F) 
C 

C CC3PUTF.S  THE  ZEEO  OECER  BESSEI  FUNCTION. 

C 

c 

C  DIfiENSION    A1  (  16)  ,A2  (11)  ,  A3  (1  1)  ,A4  (15)  ,A5(  11)  ,A6  (1  1) 

EEAL*3     A1  (16)  ,  A/  (11)  .A3  (1  1)  ,  A  U  (15)  ,A5  (1  1)  ,Ao  (1  1) 

DATA  A  1/-1.0D-  15,4.  1t-14,-1. 9  44D-12,  7.6  4S7D-11 ,-2. o792  57D-9,7.60  81 
1  6  36D-8,-  1.761S16BC8D-6.3.  24  60  2 26921 D  -5.  -4.  t>0  62o  1  bo2  06D-4  .  4  .  61  5  1800 
26946 8D-3.-. 03489376541140 5.. 156 06 710233 2097, -.37009499 38 7265, .2651 
37  8  6 132 0  3  337. -.CC 67234423 5  28 52,-3  1545  594  2  94 57 8/ 

DATA  A  2/  1.0D-  15.-5. 0 C- 1 5 . 4. 3D- 14,-4.  3D- 13, 5.  1 o 3D- 1 2 , -7 . 36 4 1 C-  1 1 ,  1. 
16  30b46D-5,-5. 1 7C5S45E-8 , 3 .0 75  1  84788 D-o , -5. 36  52  20  4  63 13D-4 ,  1 .  9  98  92Co 
29869504/ 

DATA  A3/-1.0D-15,4.0D-15,-3.3r-14,3.01D-13,-3.207D-12,4.220  1D-11,- 
17. 2719  2D- 10,  1.75/245 7D- 8, -7. 4 144 9841 D- 7,6. §385 199426D- 5, -.031  11  170 
29210674/ 

DATA  A4/1.1D- 14.-5.7  8D- 13,2. 5281D-11  ,-9 . 4242 1D-1 0 , 2 . 94 9707 D-8  ,-7  .  6 
11753781D-7, 1 . 56 E7C 1 924D- 5 ,- -  0002604 44385349,. 0  03 240270 132684, -.C29 
2175  524  806  154,  . 177709 1172 3972 6, -.66 14 43934 134 543, 1.2  37994  096  65768,- 
31. 191801 16054122, 1.2  567 17  54  12  1C53/ 

DATA  A 5/-  1.0 D-  15,  5.  OD-15,-4. 7 £-14, 4. 7D- 13 . -5. 705D- 1 2,6 . 81 6 9D- 1  1,-1 
1.8718910-9,6.  177  6  34D-8,-3.9  87  2£43D-6,. 00089898 93330 So, 2.  00 1 60  60  8 17 
22003/ 

DATA    A6/1.0D-15,-5.0D-15,3. 6D-14.-3.  26D- 1 3,  3  .  5 1 5D- 1 2, -  4. 6 86  4D- 1  1,8 
1.22919D-1C,-2.C5  557£  1D-8,  5.  1  36  fc  1526D-7  ,-9  .  6277235492D-5  , .  09  3555  574 
2139071/ 
C ZEBC    OEDEE    EESSII    JUNCTION 

N  =  0 
C 

IF     (N.  EQ.  1)     GC    TO    50 

IP    JABS  (X).GT.S. 0)     GC    TO     20 

Y=.0625*X*X-2 

E  =  0.0 

BF1=0. 0 

DC     10    1=1,16       . 

EF2=BP1 

EF1  =  B 
10  B=Y*BP1-EF2+A1  (I) 

F  =  .5*  (B-EE2) 

BET DEN 
20  Y  =  256/ (X*X) -2 

AE=ABS  (X) 

B  =  0.0 

BP1=0. 0 

DC    30     1=1,11 

BP2=BP1 

EP1=B 
30  B=Y*BP1-EF2+A2  (I) 

P=.5*  (B-EE2) 

B  =  C.O 

BF1=0. 0 

DC    40    1=1,11 

EF2=BP1 

BF1=B 
40  B=Y*BP 1-EE2+A3  (I) 

Q  =  4*  (B-BF2)/AE 

Y=AB-.  785398U23S7446 

F=.7978845608C2865*  (F*COS  (Y) -C*SIN  (Y  )  ) /SQRT  {  AB) 

EETOEN 
50  IF     (ABS  (X).GT.  8.0)    GC    TO    70 

Y=.0625*X*X-2 

B  =  C.O 

BP1=0. 0 

DC    60    1=1,15 

BF2=BP1 

BF1  =  B 
60  E  =  Y*BP1-EJ2  +  A4  (I) 

F=.0625*  (E-BP2)*X 

EETURN 
70  Y  =  256/ (X*X)-2 

AE=ABS  (X) 

B  =  0.0 

EP1=0. 0 

DC    80    1=1,11 

EF2=BP1 
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80 


90 


£P1  =  B 
B=I*BP 
P  =  -5*  ( 
B  =  0.0 
BP1=0. 
DC    90 
EE2=BP 
BP1  =  B 
E=Y*BP 
Q=4*  (B 
B=SIGN 
i=AB-2 
F=.797 
E  =  E*F 
EETUHN 
ENE 


1-EE2+A5  (I) 
E-EE2) 


0 

1=1 

1 


11 


1-EE2+A6  (I) 
-3E2)/AE 
[1.0, X) 

.35c19ii4SC1S23  5 

88U560  3C2£cS*(E*COS(Y)-C*SIN(Y))/SCEl(AB) 
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c 
c 

c™ 

c 


SUEBOUTIKI    FIE1I 


(E,  DE) 


COMPUTES    EXIT    APEFTUFE    FIELD    E1SIEIB U1ICH . 
AXI-SYflETBIC. 

CCBHON    /ABBPLE/    I II,  EIAOBS,  BE  At!  SZ  ,  00 

OR  =  0. 

D=2.*B 

IF     (D.GT.EIA.  CE.D.ll.EIAOES)      FETUBN 

EBEAD=EIABSZ/Z. 

UR=U0*EXE  (- (E/EBEAE)  **2) 

BEIUBN 

ENE 


NOTE    IT    BUST    3E 
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10 


SUEEOUIIKI    DOCS']     (E1,N) 

CCHPUTES     NCRMAIIZATICN    COKSTAM    GIVING     ETOTAL    ENERGY    IN    SEAM 


COB 
PI  = 

bc= 

ECS 
BI  = 
BIS 
DX  = 
CST 
X  =  D 
□  0  = 
SDB 
DC 
BCO 
If 
CIlL 

sua 
x=x 

TEB 
00  = 
BET 
END 


HON 
2.* 
DIA 
Q=fi 
DIA 
Q  =  R 
BOS 
=  PI 
X/2 

=6. 

10 
TX= 

=  SU 
+  DX 
P  =  S 

SQR 
UEN 


I2A,:iAOBS, 3EA^SZ,U0 


/1RBFLI. 
ARS1N  (  1.; 

i2- 

0**2 

OBS/2. 

1**2 

Q/N 
*DI 


1=1.  N 

SQET  (X) 

OTX.GT. EC.CE.ECOTX.  IT.RI)     GO    TC     10 

IEID     (HCCTX,US) 

M+0E**2 

UH*CST 

T  (El/TEBE) 
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c 

c— 

c 


SOEBOUTIhE    HTIATK     (X  ,D, HHC, AD  AT , f ) 

■-    CCMPDTE5     ATMOSfEEEIC    H7F     FUNCTION 

XE=X/D 

A  1=1.-ADAE 

DEHO=D/EHC 

F=EXP (-XE**  (1.66  67)  *  (1. -A  1*XD**  (.  3333)  )  *D£HO**  (1.6667)  ) 

RETURN 

END 
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c 

c— 

c 


SUEB03TIM    KTIoIl     ( 1  ,SIG  JIT  ,  W  i  VE,  F) 

--    CCHPDTES    JITTEE    E1F    FUNCTION 

TFI=4. *AKSIN (1 .) 
CK=TPI/H17E 

F=EXP (- (CK*X*£IGJIT) **2/8.) 

BE1UBN 

ENE 
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c 

c— - 
c— 

c 


S0E8OUTIKE    3TFEC      U,F) 

CCHPOTES    SAHDCfl    FBASE    MTF. 
FOBCTION  . 


ASSUMES     A    GAOSSIAN     CORRELATION 


CCMMON    /EC/    TEIIMT,  VARBQ,  SAVEECSCALBQ 

IEMP=1  .-EXP(-  (X/SCALEQ)  **2) 

F=EXP  (-VAF30.MIMF) 

RETURN 

ENE 
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10 


S0ERO0TIKE    BQIHI1     ( A  ,  N4  ,  DX,  E  El) 

DlfTENSIOK    A(N4) 

SD«=0. 

X=DX/2. 

DC     10    1=1, NU 

Clll    MTFEC.     (X,F) 

SDH=SUM+A  (I)  *f  *X 

X=X+DX 

CONTINUE 

BZ1=SUB*CX*6. 2621852 

EETURN 

END 
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c 

c— 

c— 

c— 

c 


10 


20 
30 


40 


SUEEOUTIM    DETWV     (A,DX,N4) 

IF    THE    BBS    WAVES    CF     PHASE    DISTCETION     ABE    N  CT    SPECIFIED, 

THEN    THIS    BOUTINE    Hill    DETERMINE    THE     APER0P2IAT2    HAVE3Q    THAT 

PBODUCES    TEZ     SEEC1FIED    TDfl.1T. 

CCHMON    /EQy    T  Eflt  I,  V  ABB  Q,  SiAV  EEC, SCALES 

DIBENSIOl>    A(N4) 

EEI0=1  ./IDFLi1T**2 

VABEQ=-AICG(HEI0) 

DVEQ=VAEEC./1C. 

IEST=. 005 

NSIGN0=1. 

CONTINUE 

CALL  BQIEEL  (  A  ,  N 4  ,  EX  ,HEL) 

D£EI=BEL-BELO 

If  (A3S  (EEEL) /EELC. LI. TEST)   GC  TO  40 

IF  (DBE1.IT.0. )  GC  TC  20 

NSIGN1=1 

GC  TO  30 

NS1GN1=-  1 

continue 

if  (nsign1.ns.  ksigno)  d  v  3q=d  vec./2- 

vabbq=vaebq+bve;*ksigni 

nsign0=n£1gn1 

GC    TO     10 

CONTINUE 

TBOPI=6. 2831852 

HAVEBQ=SCET(7AEBC)/T«OPI 

BETUBN 

END 


119 


c 
c— 

c 


10 


SUSROU1IKI    JI1      (A,DX,N4  ,  S  IG  JI1  ,  W  A72) 

CCHPUT25    AEBAY     FCE    JITTEE    .IT?. 

DIBENSIOh    A(NU) 

X=EX/2. 

DC     10    1=1,  N4 

CALL    MTFJIT     (  X  ,SIG  JI  T,  B  A  V  E,  F) 

A(I)=F 

X=X+DX 

EE1URH 

ENE 
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c 
c— 

c 


10 


SOEBOUTINE    B3      |A,IX,N4) 

— CCHPOTES     4BEA3     POE    EEAM     QCALI1S    MTF. 

DIMENSION    A (N4) 

X=DX/2. 

DO     10    1=1, HI 

CAIL    HIFEQ     (X,E) 

A(I)=F 

X=X+DX 

EETOEN 

END 
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SUEROUTINE    CH2H     ( t E I GHT , CN2 ) 

CACULATE    ATnC££fcEEIC    VERTICAL    1UE30LENCE 

--    HUFNAGELS    LATEST    SCD2L-GCCD    OMY    ABOVE    3    KILOH2TSES 

CH2=2.2*  (1.E-5*10.**  (-.3)  *H EIGHT)  **1 0. *EXP  (-HEIGHT/ 1000.) 

CH2=CH2+1-S-16*EXE (-EEIGHT/15C0.) 

CH2=CN2*2-7 


c 

C THIS  MODIFICATION  IS  OSED  to  include  turbulence  at  lower  alts 

c 

CN2=CM2+  1-2-1  4*Ei£  (- tEI  3  3  T/1  0  0  .) 

BETURN 

ESD 


122 


c 

C— 

c 


10 


SUBROUTINE    AV      |N  ,  10  ,  EHTDCI,  HT  ,  CCSW  ND  ,cosox  j,  E) 

COMPUTES    TEAT    fABI    OF    THEBEAL    E1CCMING    TEAT    CHANGES    WITH    SLANT    FATh 

COBWON    /AT.10/     EATfiC 

COSSQ=CCSSND 

DE1H= (HATeO-HT)/N 

RG=DELH/2. 

HEIGHT=HT*DELE/2. 

V0COS=V0*COSSC 

S0B=0. 

SUHAS=0. 

DC    10    1=1, H 

CAIL    ALEA     (ALA, HEIGHT) 

CALL    ALES     |ALS. HEIGHT) 

SDKAS=S0HA3+AIS+AIA 

ALCSS=EXE  (-S0S1S*EEIE/C0SCMG) 

SUK=SUM+AICSS*ALA/  ( V CCO  S  +  HG*? BICCT) 

EG=BG+EEIE 

HEIGHT=HE1GHT+C£LH 

CONTINUE 

SU«  =  Sat1*LELH 

E=SUM**2 

RETURN 

EKC 
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10 
20 


SUBROUTINE    PRTFBS     (  E  I  ,K1H  ,  W  A  V  I  ,PH) 

■     COMPOTES    THAT     EAET    CF    THERHAi    BLOOMING    PHASE    THAT    DOES    NOT     DEPEND 
•     CN    SLANT    EANGI 


COMMON 

CIMENS 

BLfl=ML 

MLHM=M 

DIASQ= 

BADSQ= 

PIU=AR 

CK=8.* 

SCALE= 

ECSQ 

EISQ 

DEIN 

XI=DEL 

YI=-DI 

1  =  11 

DC    20 

I  =  YI 

SUB=0. 

XSQ=1* 

DC    10 

ESC=Y* 

E=SCRT 

FACTOR 

IF     IBS 

FACTOR 

SD«=SO 

PH  |I,  J 

Y=Y+D2 

CGNTIN 

X=X+DE 

CGNTIN 

EE10BN 

END 


ION 
-1 

LH- 
DIA 
DIA 
SIS 
PI« 
CK* 
DIA 
DIA 
IA^ 
Vi 
A/* 

1=1 


EBFLE/    IIA,DIAOES, BEAffSZ,00 
PH  (M1K,ML) 

1 

**2 

SC/U. 

/  K  A  V I 

(2.72I-U*.4)/(1.01E5*1.4) 
/2.  )!*2 
CBS/2.) **2 

MLM 

-+DE1V2. 

,MLHB 


1  ,MLM 
XSQ 

SO 


Il.RCSC.AND. ESQ.  GT.  RISC)     CALL    FIELD     (1-  ,  FACTOR) 
ACTOr**2 


LN 
UE 
LN 
OE 


FACTOE 

SGM*SC AIE*DEIN 
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SUEROUTINE    RESEHS     (H1,J1LH,PH) 
C 

C CALCULATES    TEE     ESSIDUAL     PHASE     AFTER     TEE    MEAN, TILT,     AND    FOCUS 

C HAVE    BEEN    SEMCVEE. 

C 

DIKENSION    PH(KIK,HL) 

COMMON    /AEBFLE/     E 1  A , EIA OES, EE A ESZ  ,U0 

HLHM=MLH- 1 

BIB=HL-1 

DELN=DIA/KLM 

XI=DELN/2. 

YI=-DIA/2.+DELN/2. 

ECSQ= (DIA/2. ) **2 

BISQ= (DIACBS/2  .) **2 
C 

C  COMPOTE    NCEEAL1ZAIION     CONSTANTS 

C 

CALL    COEFS     ( A C , 1 1 , A2  ,  A3 , A  4, H L H , MLHM , DELN) 
C 

C COMPUTE    EXPANSION    COEFFICIENTS 

C 

SUK1=0. 

SDH2=0. 

SDH3=0. 

SUM4=0. 

X  =  XI 

DO    20    I=1,flLHM 

XSQ=X**2 

Y  =  YI 

EC    10    J=1,l!Lfl 

Y£Q=Y**2 

BSC=XSC+YSQ 

IF     (ESQ. GT.BCSC-OE. ESQ. LT. RISC)     GO    TC     10 

CALL    FIEIE     (SCET  (FSCJ  ,UP.) 
C 

C DETERMINE    EXPAKSICN     COEFFICIENTS    RELATIVE    TO    A    UNIFORM 

C WEIGHTING    FUNCTICK. 

C 

FI=0. 

IF    JABS  (EE). GT. 1.1-2)    FT=1. 

FACTOR=PE  (I, J)  *FT 

SUM1=SUM 1+FACTCE*X 

SDH2=SUn2+FACTCF*I 

SDM3=SUM3+PACTCE*ESC 

SUM4=SUH4+FACTCE 
10            Y=Y+D£LN 
20             X  —  X  +  DEXiN 
C ADJUST    FOE     HALF    PLANE    INTEGRATION 

S0M1=0. 

SUM2=SUM2*2. 

SUB3=SDM3»2. 

SUfl4=SUM4*2. 
C 

DELNSQ=DE1N**2 

?HMEAN  =  AC*SUM4*EELNSC. 

PHILTX=A1*SU3  1*D£LNSC 

PBTLTY=A2*SUM2*EE1NSC 

PBFOC=  (A:-*SOM3-»AU*SUMU)  *DE1NS£. 
C 

c 

C SUBTRACT    THE     M  I  A  N  ,IIIT,  ANE     FOCUS    CURVATURES. 

C 

X  =  XI 

DC    40    I=1,MLHM 

XSC=X**2 

Y  =  YI 

DO    30    J=1,MLM 

YSQ=Y**2 

BSC=XSC+YSQ 


IF    jRSQ. GT. ROSC.OE.ESQ.LT. RISC)     GO    TC    30 
PH(I,J)=FH(I,J)-FETLTX*A1  *X  -P  ETLT  Y*  A  2*  Y  -FH] 


1FOC*  (A3*ESQ4-A4)  -PH3EAN*A 

10       ••-•-    • 

30  Y=Y+DELN 

40  X=X+DELN 

RETURN 

END 
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SUBROUTINE    CCEFS     (AO ,41 , A  2, A  3 , A4 ,MLM ,HIB E,D£LN) 
C 

C COMPUTES    TILT     A8D    FOCUS     HCRMA1LZTICN     CONSTANTS 

C 

COMMON    /1EBFLI/    EI2, EIAOES, BEAKSZ ,00 

BOSQ= (DIA/2. ) **2 

RISQ= (CI5CBS/2.) **2 

DELN5^=DELN**2 

YI=-DLA/2.+DELS/2. 

XI=DELN/2. 

2  =  XI 

SUM1=0. 

SDK2=0. 

SUK3=0. 

SDMU=0. 

DO    20    I=1,MLHM 

XSC=X**2 

X=fl 

DO     10     J=1,MLM 

YSC=Y**2 

ESC=XS£+YSQ 

BSCSQ=RSC**2 

IP     (RSQ.  CT.RCSC.CE-ESQ.LT.EISC.)     GO    TO    10 

CALL    FIELD     (SC.  FT  (ESQ)  ,DH) 
C 

C DETERMINE    NORMALIZATION     COEFFICIENTS     RELATIVE    TO     A    UNIFORM 

C WEIGHTING    FUNCTION. 

C 

PT=0. 

IF     (ABS  (CE)  .GT.  1.E-2)    FT=1. 

FACTOR=FT 

SUM1=SUM 1+FACTCR 

SUM2=SUM2*FACTCF*XSQ 

SUM3=SUM3+FACTCE*YSC 

SUMU=SUM«+FACTCB*ESCSQ 
10            Y=Y+DELN 
20            X=X+DELN 
C ADJUST    FOB     EALF    PLANE    INTEGRATION 

SUM1=SUM  1*2. 

S0M2=SUM2*2. 

SUM3=SUM3*2. 

S0MU=SUM«*2. 
C 

A0  =  SQRT  (1-/ (SUE1*EELKSQ 

A1=SQRT     1./(S0M2*EELNSQ 

A2=SQBT  (1./(SUM3*DEIKSQ 

E=-SUM  VJ2UM2  +  EUK3) 

A3=SQRT (E*E/ ( (  E* S  +  5UMU- S UM 1 )  *EELNS^) 

AU=A3/B 
C 

BETUBN 

EHD 
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10 

20 

30 
C  — 


SUEROUTIKE    PHViB      |B  I  ,  MLH  ,  ETOT  11 ,  PH,  S  IGS  £0) 

•COMPUTES    THE     IfcEEKAI    BLOOKING     EHASE     VARIANCE 

•VARIANCE    IS    CCtfUTED    LIKE    THE    STREHL     RATIO    ACCORDING    TO    THE 
FIELD    AS    A    BE1GHI1NG    FUNCTION. 

COMMON    /ARBFLD/     E I  J, EIAO E £, BE ABSZ , UO 

DIMENSION    PH(Plfi,P.I) 

RCSQ=(DIA/2.)  **2 

RISQ=JDIAC3S/2.) **2 

MLB=ML-1 

MIHM=KLH-1 

DEIN=DIA/MLH 

XI=DELN/2. 

YI=-DIA/2.+XI 

X  =  XI 

SUM1=0. 

SDM2=0. 

SUM3=0. 

DC    30    I=1,MLHK 

Y=1I 

xsc=x*x 

DO    20    J=1.HLM 

ESC=XSQ+Y*Y 

IF     (ESQ-GT. RCSC.CE.ESQ.LT. RISC)     GC    T C     10 

CAIL    FIELD     (SCET  (ESQ) ,UR) 

SUM1  =  5UM  1  +  Ufi 

SUa2=SUB2+UR*EE  (I,J)  **2 

SDM3  =  SOM3  +  US*Efc  (I, J) 

CONTINUE 

Y=Y+DELN 

CONTINUE 

X=X+D2LN 

CONTINUE 

ADJUST    FOB     HALF     EIANE     INTEGRATION 

SUM1  =  SUM  1  +  2. *EEIN«*2 
SUM2=SUM2+2. *CEIN**2 
SDM3=SUMi+2. *DIIN**2 

NORMAII2E    HITH    RESPECT    TO    INTEGRA  I    OF    THE    FIELD. 

SIGSQ0  =  SUM2/SDM1-  (SUM3/SOB1)  **2 

RETURN 
END 
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c 
c 

c— 

c 


SOEROUTINE    ISC1ES     (N  ,  HI , HG, K ,  ISCAKG) 

COMPOTE    ZEXITK     ISCPLAiiATIC    ANG1E    BASED    DPON    D. 
•DEE  INI  TICK 


FHIEDS 


10 


COflHON    /A1HO/     EAlffC 

BEAL    ISCAKG 

TPI=4. *AESIN  (1 .) 

CK=TPI/H 

CKSQ=CK*CK 

EEIB= (HA1BO-aT) /H 

HEIGHT=HG+DELo/2. 

DH3=DELH/2. 

SUM=0. 

DC    10    I=1,N 

CALL    CN2E     (EEIGH1CCN2) 

S0I«=SUH  +  CN2*DHT«*1.  6(666667 

HEIGHT=KE1GHT+EEIE 

EHT=DHT+CELH 

CCKTINOE 

CNI=SOM*rELH 

EEN=2. 91 *CKS;*CNI 

ISOANG=  (6.88/EEN) **. 6*. 314 

EETDHN 

EKE 
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10 


20 
30 


40 


50 
60 

70 
80 


90 

100 

110 


120 
C 
C 
C 
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C 

C 


SUBROUTINE    ISCJIA     ( 3C  ,XNU  ,H  ,  S  I  CC  HG,  H  AVE  ,  XINT,  HT,  HG,  TISC  ,  N  4  ,  II) 

■THIS    SUBROUTINE    CCMPUTES     1HE     ISCFLANAIIC    BTF    FOR     A     FULL    AO 
-COMPENSATED    SISTER    DEVELOPED    EH    EL    PEIEI 


DI  KENS  ION  H(23),C(1S),C2(19)  ,  C4  (19) 

DIMENSION  TISC (N4) 

DIMENSION  S  f  19}  ,I«U  (20)  ,CN2  (22) 

COMMON  /ATMO/  EAISC 

IF     fJJ.GE.1l      GC    TC    3C 

NE  =  20 

N  =  1 8 

DELH=(HATBO-ET)/NE 

HEIGHT=HG*DELE/2. 

DC    10    K=1,NH 

CALL    CN2E     (  HE  I  GUI ,  CN  SQ) 

CK2  (K)  =CNSQ 

H  (K)  =HEIGHT 

HIIGHT=HEIGH1+DEIB 

PI=3.  141S926 

A13=1./3. 

A56=5./6. 

A572=5./72. 

A  143=1  4./3. 

A9127=91 ./27. 

A73=7./3. 

A53=5./3. 

A12=1./2. 

DEBI=PI/2./S 

PHI  =  DPEI/2. 

DC    20    J=1,N 

C  (J)=CCS  (PHI) 

PHI=PHI+EEBI 


)*C2  (J) 


C2|jj=cia^cij,. 

comtinu: 
JJ=JJ+  1 
xa=xnu*: 

DC    40    J; 
S<  J)  =0. 
DC    11C 


c 

CC 

JJ=JJ+  1 

XA=XN0*SECOHG/£ 
DC    40    J  =  1,N 


10     K=1,NH 


X=H  (K)  *XA 
XX=CN2  (K) 


IF     (X.LT.0.1)     GC    TO    60 

IF 

X2=X*X 


(X.LT.O.  1) 
(X. GT.  10.) 
=  X*X 


GC    tc    80 


DO    50    J=1,N 

F=1.*X*»A53-A12*  (1.+2.*X*C(J)+X2)  **  A  56- A  12*  (1.-2.*X*C(J)+X2)**A5£ 

S  (J)=S  (J)  +  F*XX 

GO    TO     IOC 

X53=X**AS3 

X2=X*X 

X4=X2*X2 

DC    70    J=1.H 

F  =  X53-A5  6*X2*  M.-A13*C2  ( J ) ) + A 572*X4*  ( 1 . -A  1  43*C2  ( J)  +A9127*C4  (J) ) 

S  (J)=S  (J) +F*XX 

GC    TO     100 

X13=X*»  (-A13) 

X73=X13/  (X*X 

DC    90    J=1.N 

F=1.-A56*X13*(1.-A13*C2(J))+A572*X7  3*(1.-A14  3*C2(J)+A9  127*C4(J)) 

SIJ)  =S  (J)  +  P*XX 

CONTINUE 

CONTINUE 

XHO  (1)  =1  ./XA 

DC    120    J=1,N 

XflO  (J+  1)  =114.  ee*DELB*S(J) 

NOW    INTEGRATE    CVEE    PHI 

CSl=B**A£3*SECCrG/UAVE**2 

SDB=0. 

DC    130    1=1, N 

SUM=SUM+EXP(-XMU  (1+1)  *CST) 
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XIHl=Sa«*DPHI*2./fI 
TISO  (II)  =XIN1 
EE1UBN 
END 
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SUEBOUTINE    FNEEfcC     ( EXD,  E  HO,  DI  t.  ,BHOU  ,  EX,  NU  ,  IHEf.TF) 
C 

C THIS    SOBFCUTINE    DETERMINES    WHA1    VALUE    CF    EHO    ROULD    PHODUCE    THE    AC 

C- CCEBEC1ED    STEEHI     AS     II    THERE     kIBE    NO     CGSPENSATION     AT    ALL 

C 

BEAL    IBBETF  (NU) 
C 

B=BHOU 

DH=BHOD/2- 

NSIGN0=-  1 

NSIGN=-1 
10  CONTINUE 

CALL    TREBEL     (IERK1F  ,  EIA  ,  N  <i,  DX  ,  SB  ,  El 

IF     (ABS  fSF-EXE)/EXD.  IT..  1)     GO    10    40 

DIF=SB-£XE 

IF     (DIF.I1-0. )     GO    TC    20 

NSIGN= 1 

GO    TO    30 
20  NSIGN=-1 

30  CONTINUE 

IF     (DIF.GE.O.)     DH=DB/2 

NSIGN0=NSIGN 

B=E-DR*N£IGN 

GC    TO    10 
40  CONTINUE 

EHO=B 

EETUBN 

END 
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c 

c— 
c— 
c 


SUEBOOTINE    TBEEEi     ( IEP.KT  F  ,DI  A  ,  NU  ,  DX  ,  IB  EI ,  E  HO) 


10 


THIS    S 
TILT    0 

BEAL    I 

PI=2.* 

SD«=0. 

EX2=DX 

X  =  EX2 

DC     10 

F1=IEE 

CALL    « 

SUK=S0 

X=X+DX 

CCNTIN 

TEEL=S 

BETUEN 

END 


UBBCOTI 
BEEECTE 

EEKTF  (N 
AESIN (1 

/2- 

1=  1,N4 
HTFJI) 
TFATK 
M+E1*F2 


HE    CCMtCTES     THE     RELATIVE     INTENSITY    DUE    TO    A    TOTA1 
I    SYETE*    FCE    TUBEULENCE    CNIY. 

-) 


X,EIA,EHO, 0. ,F2) 
*1 


UE 
0M*2.*PI*EX 
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SUBBOOTIM    TL1ISC     (D  ,OMX  ,QH  Y  ,  BISILT,  SECC  3G  ,H  TE  AKS  ,  K  ABVGD) 
C 

C THIS    SOBHOUTIKI    CCSFCTES    THE    EESIDUAI    1IIT    DOS    TC    ISOP LAN A TISK 

C 

COMMON    /AIMO/     fcATMO 

N=100 

E56  =  5./6- 

E53=5./3. 

DS=  (HATMC-HTBA5S)  /N 

H=BA3VGD+ES/2 

SUB=0. 

DC    10    1=1, N 

CALL    CN2E     (H,CK2) 

2=SECOBG*H 


ARG1= (  (D+2*0!U)**2  +  (2*0 HY) **2) **E56 

ABG2=  (  [D-Z*0MX1**2+     Z*OMY)**2)**E56 

ABG3=  fZ**2*{OHJ**2+OHY**2})  **E56 

ABGU=D**E53 

AEG5=  (  (Z*C.1X)  **2+  (D-Z*0MY)**2)  **E56 

ARG6= {  (Z*CMX)  **2+  (D  +  2*OMY)**2)**E56 
C 

SUM=SaH+  (i).*ABG4  +  £J.*ABG3-ABG2-A£G5-ABG1-ABG6)  *CN2 

H=H+DS 
10  CONTINUE 

C 

C THIS    IS    TEE    RESULTANT    2-SIGMA-E    TILT    DDE    TO    ISOPLAN ATI SM. 

C 

BESTLT=SDM*DS*SICC.1G*2.  9  1/D**2*2. 
C 

BETDBN 

ERD 
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SOEBOUTIKE    BZITCT     (1,TR,TH) 
C 

C if    HAVE    EOUNE     TEE    DETAILED     M I  f 

C HHICH     LIE    BEIS.EEK     AN    ESS    rIODII 

C HDLTIPLICATIVE     AEPRCACH.        THIS 

C AT    BETTEE    MATCHING    THE     EE5DLTS 

C AVERAGING    THE    ESS    AND    THE    « DIIIPLIC A TI VE 

C 

c 

T=JTB+TH)/2. 

EETOBN 
END 


CODE  TC  GIVE  HE5ULTS 
EGK  ELCCMING  AND  A 
S03EC0TINE  IS  AH  ATT: 
CF  TEE  KTF  CODE  BY 
BESOLTS. 
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c 

c- 
c — 
c — 
c — 
c — 


SDBROUTIRI    BLCCf     (D,£,S,T) 

— THESE    BLOCKING    IHE1    HCDELS    ABE    CURVE    FITTED    TO    3UTSHTF 
--BUNS    FOR     SIGSQ'S     GREATER     TBAN     AECUT    1.2.     THE    APERTURE 
— DISTBIBOTICNS     USEE    HAE    OBSCURATIONS    OF    . 1    X    THE    OUTER 
— EIAHETEfi.       IF    EEITIE    IS    RECUIRFE,    RECCftflEND    BASELINING 
-TO    THE    GUTSHTF    CCCI    AGAIN. 


SCflTS=SQ£T  fS) 

TEG=-. 08705+5C.ETS*2.S1485  +  .172  2*S 

THU=U2877-SQB3S*2.£U91+4.096C5*S 
TBU=1./TEU 
T=  (1.-  (D/E>  ** 
IF     (D. 

RETURN 
END 


/E)  **2>  *TEU  +  (D/B)  **2*TFG 
.£)     T=ir" 


[RG 
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SOEROOTINE    SCIKT      (SI  GXZ  ,  H  ,HTfi  A  hS  ,  HGR  ND  ,  «  AV  E) 
C 

C THIS    BOUTINE    CCMEOTES    THE    VARIANCE    OF    TEE    LOG    A3PLITUD! 

C FOB    A    ZENITH     ANGLE.       THE    EESUI1S    WILI    ONLY    BE    USED 

C IN    THE    EVENT    TEAT    FU1L    AO    IS     DTILIZEE    IN    THE    HON. 

c 

COHMON    /ATMO/     EATBC 

CK=6.28/SAVE 

HTOTAL=HAlflO-HTEASS 

DB=HTOTAL/N 

H=HTBANS+EH/2. 

C56=5./6. 

SOfl=0. 

DC    10    I=1,N 

H56=H**C56 

CALL    CN2E     (H.CN2) 

SUa=SOM+CN2*H56 

H=b+DH 
10  CCNTINOE 

CK76=CK**  (7./6.)  *EH*0-56 

SIGXZ=SUe*CK76 

BETOBH 

END 
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c 

C  — - 
c 

c— 

c 


SUEROUTINE    SIBILS     (T AMP , S IGXZ  ,  5ECOMG ) 

COMPUTE    TEE    LOG    AMPLITUDE    SCI  HI I1LAT ION    VARIANCE 

FOB    OFF-ZENITH    CONDITIONS    AND    THE    RELATIVE    IHRAJIAMC; 

REDUCTION    WHEN     EUIL    AO    IS    USED. 

TAHP=EXP  (-SIGXZ*SECCMG**  (11./C.)) 

RETURN 

END 
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